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Abstract 

We prove the complete monotonicity on (0, oo) n for suitable inverse powers 
of the spanning-tree polynomials of graphs and, more generally, of the basis 
generating polynomials of certain classes of matroids. This generalizes a result 
of Szegd and answers, among other things, a long-standing question of Lewy 
and Askey concerning the positivity of Taylor coefficients for certain rational 
functions. Our proofs are based on two ah initio methods for proving that P~P 
is completely monotone on a convex cone C: the determinantal method and the 
quadratic- form method. These methods are closely connected with harmonic 
analysis on Euclidean Jordan algebras (or equivalently on symmetric cones). 
We furthermore have a variety of constructions that, given such polynomials, 
can create other ones with the same property: among these are algebraic ana- 
logues of the matroid operations of deletion, contraction, direct sum, parallel 
connection, series connection and 2-sum. The complete monotonicity of P~@ 
for some /3 > can be viewed as a strong quantitative version of the half-plane 
property (Hurwitz stability) for P, and is also related to the Rayleigh property 
for matroids. 
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1 Introduction 



If P is a univariate or multivariate polynomial with real coefficients and strictly 
positive constant term, and /3 is a positive real number, it is sometimes of interest to 
know whether P~^ has all nonnegative (or even strictly positive) Taylor coefficients. 
A problem of this type arose in the late 1920s in Friedrichs and Lewy's study of the 
discretized time-dependent wave equation in two space dimensions: they needed the 
answer for the case P{y 1: y 2 ,y 3 ) = (1 - ^(l - y 2 ) + (1 - yi)(l - y 3 ) + (1 - y 2 )(l - y 3 ) 
and (3 = 1. Lewy contacted Gabor Szego, who proceeded to solve a generalization of 
this problem: Szego [135j showed that for any n > 1, the polynomial 

n 

P n ( yi ,...,y n ) = Y^Jlfl-yj) (1.1) 

i=i j^i 

has the property that P~^ has nonnegative Taylor coefficients for all (3 > 1/2. (The 
cases n = 1, 2 are of course trivial; the interesting problem is for n > 3.) Szego's proof 
was surprisingly indirect, and exploited the Gegenbauer-Sonine addition theorem for 
Bessel functions together with Weber's first exponential integral^] Shortly thereafter, 
Kaluza [SB] provided an elementary (albeit rather intricate) proof, but only for n = 3 
and (3 = 1. In the early 1970s, Askey and Gasper [12J gave a partially alternate 
proof, using Jacobi polynomials in place of Bessel functions. Finally, Straub [134j has 
very recently produced simple and elegant proofs for the cases n = 3,4 and (3 = 1, 
based on applying a positivity-preserving operator to another rational function whose 
Taylor coefficients are known (by a different elementary argument) to be nonnegative 
(indeed strictly positive). 

Askey and Gasper, in discussing both Szego's problem and a related unsolved 
problem of Lewy and Askey, expressed the hope that "there should be a combinatorial 
interpretation of these results" and observed that "this might suggest new methods" 
[j~2| p. 340]. The purpose of the present paper is to provide such a combinatorial 
interpretation, together with new and elementary (but we think powerful) methods 
of proof. As a consequence we are able to prove a far-reaching generalization of Szego's 
original result, which includes as a special case an affirmative solution to the problem 
of Lewy and Askey. Indeed, we give two different proofs for the Lewy-Askey problem, 
based on viewing it as a member of two different families of generalizations of the 
n = 3 Szego problem. Our methods turn out to be closely connected with harmonic 
analysis on Euclidean Jordan algebras (or equivalently on symmetric cones) |56j . 

1 These formulae for Bessel functions can be found in [1441 p. 367, eq. 11.41(17)] and [1441 
p. 394, eq. 13.3(4)], respectively. For the special case n = 3, Szego also gave a version of the proof 
using Sonine's integral for the product of three Bessel functions [1441 P- 411, eq. 13.46(3)]. Szego 
commented in his introduction [1351 p. 674] that "Die angcwcndctcn Hilfsmittel stehen allerdings in 
keinem Verhaltnis zu der Einfachheit des Satzes." ("The tools used are, however, disproportionate 
to the simplicity of the result.") 

Szego in fact proved the strict positivity of the Taylor coefficients for all n when (3 = 1, and for 
n > 4/3/(2(3 — 1) when f3 > 1/2 [135| Satze I— III] . In this paper we shall concentrate on nonnegativity 
and shall not worry about whether strict positivity holds or not. But see Remark 1 after Theorem l2.2l 
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From a combinatorial point of view, one can see that Szego's polynomial ( 11. ip is 
simply the spanning-tree generating polynomial 7b(x) for the n-cycle G = C n , 



T Cn {x 1 ,...,x n ) = \\ x i > ( L2 ) 

after the change of variables Xi = 1 — ^. This suggests that an analogous result 
might hold for the spanning-tree polynomials of some wider class of graphsj^ This 
conjecture is indeed true, as we shall show. Moreover (and this will turn out to be 
quite important in what follows), the change of variables X{ = l—yi can be generalized 
to Xi = q — yi for constants q > that are not necessarily equal. We shall prove: 

Theorem 1.1 Let G = {V,E) be a connected series-parallel graph, and let Tq(x) be 
its spanning-tree polynomial in the variables x = {x e } e( zE- Then, for all (3 > 1/2 and 
all choices of strictly positive constants c = {c e } e ^E, the function Tq(c — y)^ 13 has 
nonnegative Taylor coefficients in the variables y. 

Conversely, if G is a connected graph and there exists (3 G (0, 1) \ {|} such that 
Tq(c — y) _/3 has nonnegative Taylor coefficients (in the variables y) for all c > 0, 
then G is series-parallel. 

The proof of the direct half of Theorem 11.11 is completely elementary (and indeed 
quite simple). The converse relies on a deep result from harmonic analysis on Eu- 
clidean Jordan algebras [HZ1IE3] [23 Chapter VII], for which, however, there now exist 
two different elementary proofs |127[ 154"] |131] . 

Let us recall that a C°° function f(x\, . . . , x n ) defined on (0, oo) n is termed com- 
pletely monotone if its partial derivatives of all orders alternate in sign, i.e. 

'-'''A >- ™ 

everywhere on (0, oo) n , for all k > and all choices of indices i\, . . . , ik- Theorem ll.il 
can then be rephrased as follows: 

Theorem 1.1' Let G = {V,E) be a connected series-parallel graph, and let T^(x) 
be its spanning-tree polynomial. Then Tq 13 is completely monotone on (0, 00)^ for all 
f3 > 1/2. 

Conversely, if G = (V,E) is a connected graph and there exists (3 G (0,1) \ {|} 
such that T^f is completely monotone on (0,00)^, then G is series-parallel. 

Allowing arbitrary constants c > thus allows the result to be formulated in 
terms of complete monotonicity, and leads to a characterization that is both necessary 
and sufficient. Szego's result (or rather, its generalization to arbitrary c) extends to 
series-parallel graphs and no farther. 



2 See (|4.1j) below for the general definition of the spanning-tree polynomial Tb(x) for a connected 
graph G pi 1130]. 
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But this is not the end of the matter: we can go far beyond series-parallel graphs 
if we relax our demands about the set of (3 for which T^ 13 is asserted to be completely 
monotone. The key here is Kirchhoff 's matrix-tree theorem pUl ETJ 11071 11031 EE1 S3 
ESI 11501 11041 [T] , which shows how spanning-tree polynomials can be written as de- 
terminants. This line of thought suggests that complete monotonicity of P _/3 might 
hold more generally for the homogeneous multiaffme polynomials arising from deter- 
minants of the type studied in [UU Section 8.1]. This too is true; in fact, such a result 
holds for a slightly more general class of polynomials that need not be multiaffine. 
We shall prove, once again by elementary methods: 

Theorem 1.2 Let A±, . . . , A n (n > 1) bemxm real or complex matrices or hermitian 
quaternionic matrices, and let us form the polynomial 

P(x 1 ,...,x n ) = detfy^XiAA (1.4) 

in the variables x = (x\, . . . ,x n ). [In the quaternionic case, det denotes the Moore 
determinant: see Appendix^] Assume further that there exists a linear combination 
of A\, . . . , A n that has rank m (so that P ^ 0). 

(a) If A\, . . . , A n are real symmetric positive-semidefinite matrices, then P _/3 is 
completely monotone on (0, oo) n for (3 = 0, |, 1, |, . . . and for all real (3 > (m — l)/2. 

(b) If A 1: ... , A n are complex hermitian positive-semidefinite matrices, then P _/3 is 
completely monotone on (0, oo) n for (3 = 0,1,2,3,... and for all real (3 > m— 1 . 

(c) If Ai, . . . , A n are quaternionic hermitian positive-semidefinite matrices, then 

is completely monotone on (0, oo) n for (3 = 0,2,4,6,... and for all real 
13 > 2m - 2. 

These curious conditions on (3 are not just an artifact of our method of proof; they 
really are best possible. They can be better understood if we take a slightly more 
general perspective, and define complete monotonicity for functions on an arbitrary 
open convex cone C in a finite-dimensional real vector space V (see Section [2]). We 
then have the following result that "explains" Theorem 11.21 

Theorem 1.3 (a) Let V be the real vector space Sym(m,IR) of real symmetric m x 
m matrices, and let C C V be the cone Il m (M) of positive- definite matrices. 
Then the map A i— > (det A)~@ is completely monotone on C if and only if (3 e 
{0, i 1, §,...} U [(m - l)/2, oo). Indeed, if (3 £ {0, §, 1, §,...} U [(m - l)/2, oo), 
then the map A y (det A) -13 is not completely monotone on any nonempty open 
convex subcone C'CC. 

(b ) Let V be the real vector space Herm(m, C) of complex hermitian mxm matrices, 
and let C C V be the cone Il m (C) of positive-definite matrices. Then the map 
A i — y (det A)^ 13 is completely monotone on C if and only if (3 G {0, 1, 2, 3, . . .} U 
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[m — 1, oo). Indeed, if (3 {0, 1, 2, 3, . . .} U [m — 1, oo), then the map A \-t 
(detA)' 13 is not completely monotone on any nonempty open convex subcone 
C C C. 

(c) Let V be the real vector space Herm(m, H) of quaternionic hermitian m x m 
matrices, and let C C V be the cone II m (EI) of positive-definite matrices. 
Then the map A i— > (detv4) _/3 is completely monotone on C if and only if 
/3 G {0,2,4,6,. . .} U [2m - 2,oo). Indeed, if /3 £ {0,2,4,6,. . .} U [2m - 2,oo), 
then the map A i-> (det A)' 13 is not completely monotone on any nonempty open 
convex subcone C C C. 

In particular, if the matrices A±, . . . ,A n together span Sym(m, M), Herm(m, C) or 
Herm(m, H) [so that the convex cone they generate has nonempty interior], then 
the determinantal polynomial (11.41) has P _/3 completely monotone on (0, oo) n if and 
only if (3 belongs to the set enumerated in Theorem ll.2|[^l 

The proof of the "if" part of Theorem 11.31 is completely elementary, but the "only 
if" part again relies on a deep result from harmonic analysis on Euclidean Jordan 
algebras, namely, the characterization of parameters for which the Riesz distribution 
is a positive measure (Theorem 16.81 below; but see [127[ 154] |131] for elementary 
proofs). In fact, when Theorem 11.31 is rephrased in this latter context it takes on a 
unified form: 

Theorem 1.4 Let V be a simple Euclidean Jordan algebra of dimension n and rankr, 
with n = r + |r(r — 1); let Q C V be the positive cone; and let A: V — > K. be the 
Jordan determinant. Then the map x i— > A(x) _/3 is completely monotone on Q if 
and only if (3 e {0, \ , . . . , (r - l)f } or (3 > (r - l)f . Indeed, if (3 £ {0, f , . . . , (r - 
l)f } U ((r — 1)|, oo), then the map x i-> A(x)~ 13 is not completely monotone on any 
nonempty open convex subcone Q' <ZQ. 

We shall see that Theorem 11.41 is essentially equivalent to the characterization of 
parameters for which the Riesz distribution is a positive measure. The set of values 
of f3 described in Theorem 11.41 is known as the Gindikin-Wallach set and arises in a 
number of contexts in representation theory (671 H23[ I143[ |95| l54"l [551 156] - 

A special case of the construction (11.41) arises |40l Section 8.1] when B is an 
m x n real or complex matrix of rank m, and we set -P(x) = det (BXB*), where 
X = diag(x!, . . . , x n ) and * denotes hermitian conjugate. Then the matrix Aj in (11.4)) 
is simply the outer product of the ith column of B with its complex conjugate, and so 
is of rank at most 1; as a consequence, the polynomial P is multiaffine (i.e., of degree 
at most 1 in each variable separately) Jf| In particular, we can obtain the spanning-tree 
polynomial Tq(jc) for any connected graph G by letting B be the directed vertex-edge 

3 Branden [28] has very recently used this latter fact to determine the exact set of a € R for 
which the a-permanent |138j is nonnegative on real symmetric (resp. complex hermitian) positive- 
semidefinite matrices. 

4 See Proposition 16.31 below. 
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incidence matrix for an arbitrarily chosen orientation of G, with one row (correspond- 
ing to an arbitrarily chosen vertex of G) deleted; then the matrix-tree theorem tells 
us that T G (x) = det(BXB T ). Applying Theorem 11.2( a). we obtain: 

Corollary 1.5 Let G = (V,E) be a connected graph with p vertices, and let 7b (x) 
be its spanning-tree polynomial. Then Tq is completely monotone on (0, oo) E for 
P = 0, §,... and for all real (3 > (p - 2)/2. 

Likewise, we can apply Theorem 11.2( b) to the elementary symmetric polynomial 

E 2 ,a{xi, £2,23,24) = 2122 + 2123 + 21X4 + 2223 + 22X4 + 2324, (1.5) 

which can be represented in the form (11.41) with 

Al = (0 0)' A2= (o l)' A3= (l })' A4= (e- /3 'T 3 ) ^ 

or equivalently as ^2,4 (x) = det(5X£>*) with 5 = f J J j e^/ 3 ! ' ^ 6 °^^ a ^ n: 

Corollary 1.6 T7ie function is completely monotone on (0, oo) 4 if and only if 
j3 = or (3 > 1 . In particular, the function 

( £ (1-2/0(1-%)) (1-7) 

\i<«<i<4 / 

/ias nonnegative Taylor coefficients for all j3 > 1 . 

Indeed, the "if" part can be corroborated by an explicit Laplace-transform formula 
for E^4 for (3 > 1: see (16.161) below. The "only if" follows from the observation made 
after Theorem 11.31 since the matrices A\, . . . , A4 in (I1.6P span Herm(2, C). 

The second sentence of Corollary 11.61 answers in the affirmative a question posed 
long ago by Lewy [121 p. 340], of which Askey remarks that it "has caused me many 
hours of frustration" [TTJ, p. 56] @ (See also the recent discussion in [85].) Indeed, 
Lewy's question concerned only j3 = 1, and made the weaker conjecture that the 



5 Askey [HI p. 56] comments that, in his view, 

So far the most powerful method of treating problems of this type is to translate them 
into another problem involving special functions and then use the results and methods 
which have been developed for the last two hundred years to solve the special function 
problem. So far I have been unable to make a reduction in [Lewy's problem] and so 
have no place to start. 

But he immediately adds, wisely, that "it is possible to solve some problems without using special 
functions, so others should not give up on [Lewy's problem]." 
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function (11.71) multiplied by (4— y\ —2/2—2/3—2/4) 1 has nonnegative Taylor coefficients. 
This latter factor is now seen to be unnecessary^ 

Similarly, Theorem 11.2( c) applied to the quaternionic determinant det^ ^ 

ab — qq for a, b G K and q G H, with Al> . . . , A4 as above and 

As = (- ?Y A R =(l q A (1.8a) 



g 5 1 J ' u V* 1 

1 V3. V6. 1 \/3. V6. VTo, n QU , 

O5 = Z 7, Og = 2 7 ft, (1.8b) 

y 2 6 3 2 6 12 J 4 ' v 7 

yields an analogous result for the elementary symmetric polynomial 

E 2fi {x 1 , . . . ,x 6 ) = xix 2 + X1X3 + . . . + x 5 x 6 (1.9) 

of degree 2 in six variables^: 

Corollary 1.7 The function E^q is completely monotone on (0, oo) 6 if and only if 
(3 = or > 2. 

Corollaries 11.51 and 11.61 are in fact special cases of a much more general result 
concerning the basis generating polynomials Sm(x) of certain classes of matroids. 
(We stress that no knowledge of matroid theory is needed to understand the main 
arguments of this paper; readers allergic to matroids, or simply unfamiliar with them, 
can skip all references to them without loss of logical continuity. Still, we think 
that the matroidal perspective is fruitful and we would like to make some modest 



6 Ismail and Tamhankar [3J>] p. 483] mistakenly asserted that "the early coefficients in the power 
series expansion of 

{(I - r)(l - s) + (1 - r)(l - t) + (1 - r)(l - u) + (1 - s)(l - t) + (1 - s)(l - u) + (1 - t)(l - u)}- 1 

are positive but the later coefficients do change sign" , arguing that this is "because Huygen's [sic] 
principle holds in three-space." Huygens' principle indeed suggests that the coefficients approach 
zero, as Askey and Gasper [121 p. 340] observed; but this in no way contradicts the nonnegativity 
of those coefficients. 

7 If we define q 3 = 1 and q 4 = e~"/ 3 [cf. (|1.6|) ]. then q 3 , q 4 , q$,q§ are quaternions satisfying 

f 1 if i = j 
Re{q ^ ] = \l/2 iii^j 

or equivalently \qi\ 2 = 1 and \qi — q 3 ] 2 = 1 for all i 7^ j. From this it easily follows that 

, , ( xi + x 3 + xi + x 5 + x 6 x 3 q 3 + x 4 q 4 + x 5 q 5 + x 6 q 6 \ „ , , 
det _ . , . . = E 2 fi(xi,. ■ ■ ,x e ) . 

\ x 3 q 3 + x 4 q 4 + x 5 q 5 + x 6 q 6 x 2 + x 3 + x 4 + x 5 + x 6 J 
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propaganda for itjj) So let M be a matroid with ground set E, and let B(M) be its 
set of bases; then the basis generating polynomial of M is, by definition, Bj^ipc) = 
x s , where x = {x e } e( zE is a family of indeterminates indexed by the elements 

SeB(M) 

of M, and we have used the shorthand x s = [] i e . 

Now let B be an arbitrary m x n real or complex matrix of rank m, and define 
-P(x) = det(-BX5*). Then, as discussed previously, Theorem 1 1.2f a or b) applies to P 
and gives a sufficient condition for P^ 13 to be completely monotone. On the other 
hand, the Cauchy-Binet formula gives 

P(x) = det(BXB*) = |detS*s| 2 a; s (1.10) 

S C [n] 
\S\ = m 

where B*s denotes the submatrix of B with columns S. Since det B*$ 7^ if and only 
if the columns S of B are linearly independent, we see that P is a weighted version 
of the basis generating polynomial for the matroid M = M[B] that is represented by 
B (this matroid has rank m). In particular, a matroid is said to be real-unimodular 
(resp. complex-unimodular) if it has a real (resp. complex) representing matrix B, 
with a number of rows equal to its rank, such that | det B*s | 2 £ {0,1} for all So 
In this case the basis generating polynomial is precisely 5m(x) = det(BXB*). We 
thereby obtain from Theorem ll.2( a.b) the following result: 

Corollary 1.8 Let M be a matroid of rank r on the ground set E, and let Bm(x) be 
its basis generating polynomial. 

(a) If M is a regular [= real-unimodular] matroid, then Bjf is completely monotone 
on (0, oo) £ for (3 = 0, |, 1, |, . . . and for all real (3 > (r — l)/2. (This holds in 
particular if M is a graphic or cographic matroid, i.e. for the spanning-tree or 
complementary- spanning-tree polynomial of a connected graph.) 

(b) If M is a complex-unimodular matroid [= sixth-root- of-unity matroid], then B^f 
is completely monotone on (0, oo) E for (3 — 0, 1, 2, 3, . . . and for all real (3 > r— 1. 

In particular, by specializing (a) to a graphic matroid M(G) we recover Corollary 1 1.5[ 
and by specializing (b) to the uniform matroid 1/2,4 we recover Corollary 11.61 

We have also proven a (very) partial converse to Corollary ll.81 which concerns the 
cases of rank-r n-element simple matroids in which the matrices Ai, . . . , A n together 
span Sym(r, R) or Herm(r, C): see Corollary 16.91 below. 

8 See [HOj for background on matroid theory, and [JD] for background on basis generating polyno- 
mials. In interpreting Corollarv ll.8l below. please note that if G = (V, E) is a graph with k connected 
components, then the graphic matroid M(G) has rank |V| — k, while the cographic matroid M*(G) 
has rank \E\ — \V\ + k. Note also that the equivalence of "complex-unimodular matroid" with 
"sixth- root-of- unity matroid" is proven in j 10. Theorem 8.9]. 

9 This is not the usual definition of real-unimodular/complex-unimodular, but it is equivalent to 
the usual definition by virtue of (40J Proposition 8.6]. 
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Remark. There is also an analogue of Corollary 11.81 in the quaternionic case. 
Recall first the quaternionic analogue of the Cauchy-Binet formula fll.lOp [Propo- 
sition |A~Hg) below]: if B is an m x n quaternionic matrix, then P(xi, . . . ,x n ) = 

det(BXB*) = det ( XiAi J is well-defined for X = diag(xi, . . . , x n ) with x±, . . . ,x n 
real and equals the polynomial Yl det [5*5 (5*5)*] x s . (Note, by contrast, that 

S C [n] 
\S\ = m 

det B*s is in general meaningless because B+s need not be hermitian.) We can then 
define a matroid M to be quaternionic-unimodular if its basis generating polynomial 
can be represented in this way, i.e. if it has a quaternionic representing matrix B, 
with a number of rows equal to its rank, such that det [B+siB+s)*] €E {0, 1} for all S. 
For such matroids M, Theorem 11.2( c) implies that B^f is completely monotone on 
(0, oo)^ for (3 = 0, 2, 4, 6, . . . and for all real (3 > 2r - 2. 

A deeper study of quaternionic-unimodular matroids would be of interest. For 
instance, is the class of quaternionic-unimodular matroids closed under duality? Or 
even under contraction? (The class is obviously closed under deletion.) Which uni- 
form matroids U TtU are quaternionic-unimodular? 

A different notion of "quaternionic-unimodular matroid" has been introduced re- 
cently by Pendavingh and van Zwam [117j . It is not clear to us what is the relation 
between their notion and ours. □ 



Of course, £2,4 and £"2,6 are quadratic forms in the variables x, as is the polynomial 
£2,3 arising in the n = 3 Szego problem. This suggests that it might be fruitful to 
study more general quadratic forms. We shall prove, by elementary methods: 

Theorem 1.9 Let V be a finite-dimensional real vector space, let B be a symmetric 
bilinear form on V having inertia (n + ,n_,no), and define the quadratic form Q(x) = 
B(x,x). Let C C V be a nonempty open convex cone with the property that Q(x) > 
for all x G C . Then n + > 1, and moreover: 

(a) If n + = 1 and n__ = 0, then is completely monotone on C for all (3 > 0. 
For all other values of (3, Q" 13 is not completely monotone on any nonempty 
open convex subcone C C C . 

(b) Ifn + = 1 andn_ > 1, then is completely monotone on C for (3 = and for 
all f3 > (n_ — l)/2. For all other values of (3, Q' 13 is not completely monotone 
on any nonempty open convex subcone C C C . 

(c) If n + > 1, then is not completely monotone on any nonempty open convex 
subcone C C C for any f3 7^ 0. 
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Theorem 11.91 follows fairly easily from the classic work of Marcel Riesz [121] (see 
also J5TJ and [56j Chapter VII]) treating the case in which B is the Lorentz form on 



I". 



B(x,y) = xiyi - x 2 y 2 - ■ ■ ■ - x n y n , (1.11) 

and C is the Lorentz cone (= forward light cone) {x G M. n : x\ > y/a%~+777+ x^}. 
We are able to give a completely elementary proof of both the sufficiency and the 
necessity; and we are able to give in case (b) an explicit Laplace-transform formula 
for (see Proposition 17.61) . 

Specializing Theorem 11.91 with V = IR n and C = (0, oo) n to the degree-2 elemen- 
tary symmetric polynomials 

E 2 ,n{xi,...,x n ) = ^ X i X i ' ( L12 ) 

l<i<j<n 

we obtain: 

Corollary 1.10 The function is completely monotone on (0, oo) ra if and only if 
ft = or ft > (n-2)/2. 

By this method we obtain an alternate proof of Corollaries 11.61 and 11.71 — hence in 
particular a second solution to the Lewy-Askey problem — as well as of Szego's |135] 
original result in the case n = 30 We also obtain an explicit Laplace-transform 
formula for (see Corollary 17. 8p . 

Corollaries ll.6[ 11.71 and 11.101 lead naturally to the following question: If we write 
E r>n for the elementary symmetric polynomial of degree r in n variables, 

Er t n{.Xlj ■ ■ ■ j 2?n) ^ Xi^X^ ' ' ' X% r (1. 13) 

l<il<i2<---<ir-<fl 

(where we set E ^ n = 1), then for which ft > is E~@ completely monotone on 
(0, oo) n ? The cases r = 0, 1 and n are trivial: we have complete monotonicity for all 
ft > 0. Our results for the cases r = n — 1 (Theorem 11.11 specialized to cycles C n ) 
and r = 2 (Corollary II. 10p . as well as numerical experiments for (r, n) = (3,5), (3,6) 
and (4,6), suggest the following conjecture: 



10 In fancy language — which is, however, completely unnecessary for understanding our proofs 
— our "determinantal" proof of Corollarv ll.6l is based on harmonic analysis on the cone of positive- 
definite m x m complex hermitian matrices specialized to m = 2, while our "quadratic form" proof 
is based on harmonic analysis on the Lorentz cone in K™ specialized to n = 4. The point here is 
that the Jordan algebra Herm(2,C) ~ M x K 3 can be viewed as a member of two different families 
of Jordan algebras: Herm(m,C) and K x M n_1 [SHI P- 98]. Likewise, our "determinantal" proof 
of Corollary 11.71 is based on harmonic analysis on the cone of positive-definite m x m quaternionic 
hermitian matrices specialized to m = 2, while our "quadratic form" proof is based on harmonic 
analysis on the Lorentz cone in 1" specialized to n = 6; and we have the isomorphism of Jordan 
algebras Herm(2, H) ~ M x ffi. 5 [56l p. 98] . And finally, our "determinantal" proof of the n = 3 Szego 
result is based on harmonic analysis on the cone of positive-definite m x m real symmetric matrices 
specialized to m — 2, while our "quadratic form" proof is based on harmonic analysis on the Lorentz 
cone in l n specialized to n = 3; and we have Sym(2, M)~lxl 2 [56j p. 98]. 
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Conjecture 1.11 Let 2 < r < n. Then E r % is completely monotone on (0, oo) n if 
and only if (3 = or (3 > (n — r)/2. 

However, we have been unable to find a proof of either the necessity or the sufficiency. 

We remark that the elementary symmetric polynomial E rjTl is the basis generat- 
ing polynomial of the uniform matroid U r>n . So Corollary 11.101 and Conjecture 11.111 
concern the same general subject as Corollary II .8[ namely, complete monotonicity for 
inverse powers of the basis generating polynomials of matroids. 



In summary, we have two ab initio methods for proving, given a polynomial P 
and a positive real number (3, that P _/3 is completely monotone on (0, oo) n [or more 
generally on a convex cone C]: 

(a) The determinantal method (Theorems 11.21 [TUl and [TT^l see Section [6]). 

(b) The quadratic-form method (Theorem 11.91 see Section [7]). 

Interestingly, these two methods can be viewed as versions of the same construction, 
involving the determinant on a Euclidean Jordan algebra and the Laplace-transform 
representation of its inverse powers [56j Chapters II- VII]. We discuss this connection 
in Sections 16.31 and 17.21 

In addition to these two ab initio methods, we have a variety of constructions that, 
given such polynomials, can create other ones with the same property (see Section [3]). 
Among these are algebraic analogues of the graph (or matroid) operations of deletion, 
contraction, direct su nS, parallel connection, series connection and 2-sum (but not 
duality). By combining these operations with our ab initio proofs, we are able to 
prove the complete monotonicity of Tq for some values of (3 beyond those covered 
by Corollary 11.51 

Proposition 1.12 Fix p > 2, and let G = {V,E) be any graph that can be obtained 
from copies of the complete graph K p by parallel connection, series connection, direct 
sum, deletion and contraction. Then Tq is completely monotone on (0, oo)^ for 
(3 = 0, |,1,|,... and for all real (3 > (p-2)/2. 

In particular, the case p = 3 covers series-parallel graphs; this is essentially our 
proof of the direct half of Theorem ll.il . We also have versions of this proposition for 
matroids: see Propositions 16 . 1 01 and 16 . lTI b elow . Finally, in Propositions I6.12l and l6. 131 
we give excluded-minors characterizations of the class of graphs / matroids handled by 
Propositions 11.121 and 16.101 respectively. 

But even in the graphic case, we are still far from having a complete answer to 
the following fundamental problem: 

11 By "direct sum" of graphs we mean either disjoint union ("0-sum") or gluing at a cut vertex 
("1-sum"). Both of these operations correspond to the direct sum of matroids. 
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Problem 1.13 Given a graph G = (V,E), for which real numbers (3 > is the 
function completely monotone on (0,oo) E ? 

This question can be rephrased usefully as follows: 

Problem 1.13' For each (3 > 0, characterize the class Gp of graphs for which Tq 
is completely monotone on (0,oo) E . 

We will show in Section 14.11 that the class Gp is closed under minors — so that it 
can be characterized by listing the excluded minors — and under parallel connection. 
Furthermore, it is closed under series connection when (but only when) (3 > 1/2. 
In this paper we have solved Problem 11.131 in a few cases: 

• For (3 G {-, 1, §, • ■ •}, Gp = all graphs. See Corollary 11.8( a). 

• For (3 G (0, \),Gp = graphs obtained from forests by parallel extension of edges 
(= graphs with no K 3 minor). See Theorem 15.11 

• For (3 G (|, 1), Gp = series-parallel graphs (= graphs with no K4 minor). See 
Theorems 11.11 and 15.31 

So the first unsolved cases are (3 G (1, |): Might it be true that Gp = all graphs with 
no K 5 minor? Or might there exist, alternatively, other excluded minors? We have 
been thus far unable to determine the complete monotonicity of Tq 13 for the cases 
G = W4 (the wheel with four spokes) and G = K 5 — e (the complete graph K 5 with 
one edge deleted). Indeed, for k > 2 we do not even know the answer to the following 
question: 

Question 1.14 Fix an integer k > 0. Must we have Gp = Gp> whenever f3,(3' G 
(k k+i\g 

V 2' 2 / " 

Let us mention, finally, an alternative approach to "converse" results that we have 
not pursued, for lack of competence. When P _/3 does not have all nonnegative Taylor 
coefficients, this fact should in most cases be provable either by explicit computation 
of low-order coefficients or by asymptotic computation of suitable families of high- 
order coefficients (computer experiments can usually suggest which families to focus 
on). This type of multivariate asymptotic calculation has been pioneered recently by 
Pemantle and collaborators |114| 1115] II 161 [1131119] an d involves some rather nontrivial 
algebraic geometry/topology. In fact, Baryshnikov and Pemantle [HI Section 4.4] 
have recently used their method to study the asymptotics of the Taylor coefficients 
of P~P for the Szego polynomial (II .ip with n — 3, but thus far only for (3 > 1/2 1^1 
It would be interesting to know whether this analysis can be extended to the case 
(3 < 1/2, thereby providing an explicit proof that some of the Taylor coefficients 
are asymptotically negative. More generally, one might try to study the elementary 
symmetric polynomials E r>n : after the n = 3 Szego case £7 2 ,3, the next simplest would 
probably be the Lewy-Askey case £7 2 ,4 [i.e., fll.7]) ] . 



12 The formula in their Theorem 4.4 has a misprint: the power —1/2 should be j3 — I. 
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Let us recall that a polynomial P with complex coefficients is said to have the half- 
plane property 001 HE [1391 123 123 IH21 123 EH] if either P = or else P(x 1} . . .^c n ) ^ 
whenever complex numbers with strictly positive real parto We 

shall show (Corollary 12.31 below) that if P is a polynomial with real coefficients that 
is strictly positive on (0, oo) n and such that P _/3 is completely monotone on (0, oo) n 
for at least one /3 > 0, then P necessarily has the half-plane property (but not 
conversely). The complete monotonicity of P~^ can therefore be thought of as a 
strong quantitative form of the half-plane property. In particular, it follows that the 
determinantal polynomials considered in Theorem 11.21 have the half-plane property 
- a fact that can easily be proven directly (Corollary 16.51 below). The same is true 
for the quadratic polynomials considered in Theorem 11.91 see [10j Theorem 5.3] and 
Theorem 17.41 below. 

Complete monotonicity is also connected with the Rayleigh property [H] for ma- 
troids and, more generally, for multiaffine polynomials. Let us say that a function / 
is completely monotone of order K if the inequalities (11.31) hold for < k < K. Thus, 
a function is completely monotone of order (resp. 1) if and only if it is nonnegative 
(resp. nonnegative and decreasing). A function is completely monotone of order 2 if, 
in addition, d 2 f /dxidxj > for all Specializing this to / = P~^ where P is a 
polynomial, we obtain 

d 2 P dP dP 

P-rr-rr- < ^ + for all U- 1.14 

OXiOXj OXi OXj 

If P is multiaffine, then d 2 P/dxf = 0, so it suffices to consider the cases i ^ j. The 
inequality (11. 14j) is then a generalization of the Rayleigh (or negative-correlation) 
inequality in which an extra constant C = (3 + 1 is inserted on the right-hand side. 
(The ordinary Rayleigh property corresponds to (3 I 0, hence to taking / = — logP 
and omitting the k = condition.) It would be interesting to know whether the 
combinatorial consequences of the Rayleigh property — such as the matroidal-support 
property |14U] - - extend to the C-Rayleigh property for arbitrary C < oo. It would 
also be interesting to extend the results of the present paper, which address complete 
monotonicity of order oo, to complete monotonicity of finite orders K. In what way 
do the conditions on become i^-dependent? 



Finally, let us remark that the literature contains some other examples of multi- 
variate polynomials P for which P - ^ has all nonnegative Taylor coefficients, for some 
specified set of numbers /3. For instance, Askey and Gasper [13J showed that this is 
the case for 

P(x,y,z) = l--(x + y + z) + -xyz (1.15) 

13 A polynomial P ^ with the half-plane property is also termed Hurwitz stable. 
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whenever f3 > (vT7 — 3)/2 ~ 0.561553; Gillis, Reznick and Zeilberger [65] later gave 
an elementary proof. Likewise, Koornwinder [92J proved this for 

P(x,y,z,u) = 1 — -(x + y + z + u) + -(xyz + xyu + xzu + yzu) — xyzu (1.16) 

whenever /3 > 1; an elementary proof later emerged from the combined work of Ismail 
and Tamhankar |86j and Gillis-Reznick-Zeilberger [65] . 

It turns out that these two examples also have a combinatorial interpretation: not 
in terms of the spanning-tree polynomial T G (x) , but rather in terms of the connected- 
spanning-subgraph polynomial [X 24 "] 1130] 

c G (v) = e n^- 

ACE e£A 
(V, A) connected 

which has T G (x) as a limiting case: 

T G (x) = limA" (|v ' | - 1) C G (Ax) . (1.18) 

A— s>0 

If we specialize to G — C n and make the change of variables = — A(l — Zi) with 
< A < n, it then turns out that the Askey-Gasper polynomial (11.151) corresponds 
to the case n = 3, A = 1, while the Koornwinder polynomial (I1.16P corresponds to 
the case n = 4, A = 2. On the other hand, in the limit A — > we recover a multiple 
of the Szego polynomial (II. ip : this is simply a special case of (11.181) . 

In the same way that the complete monotonicity of T^ 13 is a strong quantitative 
form of the half-plane property, it turns out that the nonnegativity of Taylor coef- 
ficients of Pq^ x in these examples is a strong quantitative form of the multivariate 
Brown-Colbourn property (or more precisely, the multivariate property BCa) dis- 
cussed in [1241 [130J. But it seems to be a difficult problem to determine the set of 
pairs (A, (3) for which Pq x has nonnegative Taylor coefficients, even in the simplest 
case G = C3. We have some partial results on this problem, but we leave these for a 
future paper. 



The plan of this paper is as follows: In Section [2] we define complete monotonic- 
ity on cones and recall the Bernstein-Hausdorff-Widder-Choquet theorem; we also 
prove a general result showing that complete monotonicity of P^ 13 on a cone C G V 
implies the nonvanishing of P in the complex tube C + iV. In Section [3] we dis- 
cuss some general constructions by which new polynomials P with P _/3 completely 
monotone can be obtained from old ones. In Section H] we apply these constructions 
to the spanning-tree polynomials of graphs and the basis generating polynomials of 
matroids. In Section [5] we analyze the series-parallel case and prove Theorem 11.11 In 
Section |6] we present the determinantal construction and prove Theorems 1 1 . 2 \ [T731 and 
11.41 In Section Owe present the quadratic- form construction and prove Theorem 11.91 
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In Section [8] we present briefly the theory of positive-definite functions (in the semi- 
group sense) on convex cones — which is a close relative of the theory of completely 
monotone functions — and its application to the class of cones treated here. In Ap- 
pendix [X] we recall the definition and main properties of the Moore determinant for 
hermitian quaternionic matrices. 

We have tried hard to make this paper comprehensible to the union (not the inter- 
section!) of combinatorialists and analysts. We apologize in advance to experts in each 
of these fields for boring them every now and then with overly detailed explanations 
of elementary facts. 

2 Complete monotonicity on cones 

In the Introduction we defined complete monotonicity for functions on (0, oo) n . 
For our later needs (see Sections [6] and [7]), it turns out to be natural to consider 
complete monotonicity on more general open convex cones CcK". This is a genuine 
generalization, because for n > 3, an open convex cone is not necessarily the affine 
image of a (possibly higher-dimensional) orthant, i.e. it need not have "flat sides": 
an example is the Lorentz cone {x £ W 1 : x\ > \/ x\ + . . . + x^} in dimension n > 3. 

Definition 2.1 Let V be a finite- dimensional real vector space, and let C be an open 
convex cone in V. Then a C°° function f: C — y R is termed completely monotone if 
for all k > 0, all choices of vectors u\, . . . , u& £ C , and all xGC, we have 



where D u denotes a directional derivative. A function f is termed conditionally com- 
pletel y m onotone if the inequality ^2.1\) holds for all k > 1 but not necessarily for 



Of course, if the inequality (12. ip holds for all u 1; . . . , in some set S, then by linearity 
and continuity it holds also for all u 1; . . . , in the closed convex cone generated by 
S. This observation also shows the equivalence of Definition I2.1[ specialized to the 
case V = R n and C = (0, oo) n , with the definition given the Introduction. 



1 The terminology "conditionally completely monotone" is new, but we think it felicitous: it is 
chosen by analogy with "conditionally positive definite matrix" [21] [18] , with which this concept is 
in fact closely related [751 US US 1201 121] ■ ( Warning: The book [H] uses the term "negative definite" 
for what we would call "conditionally negative definite".) 

Please note that if / is bounded below, then / is conditionally completely monotone if and only if 
there exists a constant c such that / + c is completely monotone. But there also exist conditionally 
completely monotone functions that are unbounded below (and hence for which such a constant c 
cannot exist): examples on (0, co) are f(x) = —(a + x) a and f(x) — — log(a + a;) with a > and 
< a < 1. 

Of course, it follows immediately from the definition that / is conditionally completely monotone 
if and only if —D u f is completely monotone for all vectors u g C. In the multidimensional case 
this seems rather difficult to work with; but in the one-dimensional case C = (0, oo) it says that / 
is conditionally completely monotone if and only if — /' is completely monotone. 





k = 
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If T: (Vi,Ci) —> (V2, C2) is a positive linear map (i.e., a linear map T: V% —> V2 
satisfying T[Ci] C C 2 ) and /: C 2 — > K. is completely monotone, then it is easily 
seen that / o T: C\ — > M. is completely monotone. Conversely, if / o T is completely 
monotone for all positive linear maps T: (M n , (0, oo) n ) — > (V2, C2) /or arbitrarily large 
n, then / is completely monotone: for if (12. ip fails for some k, then we can take n = k 
and Te.i = Uj (where e, is the ith coordinate unit vector in M. n ) and foT will fail one 
of the /cth-order complete-monotonicity inequalities. 

Let us next recall some elementary facts. If / is completely monotone, then 
/(0 + ) = lim a! _ > .o,a:eC / (x) exists and equals sxvp xeC f(x), but it might be +00. The 
product of two completely monotone functions is completely monotone. If / is com- 
pletely monotone and $: [0, 00) — > [0, 00) is absolutely monotone (i.e. its derivatives 
of all orders are everywhere nonnegative) , then $ o / is completely monotone. If / 
is conditionally completely monotone and $: (— 00, 00) — > [0, 00) is absolutely mono- 
tone, then $ o / is completely monotone. (In particular, this occurs when $ is the 
exponential function.) Finally, a locally uniform limit of a sequence of completely 
monotone functions is completely monotone. 

The fundamental fact in the theory of completely monotone functions on (0, 00) 
is the Bernstein-Hausdorff-Widder theorem |146] : A function / defined on (0, 00) is 
completely monotone if and only if it can be written in the form 

POO 

f(x) = / e~ tx dfi(t) (2.2) 

where fj, is a nonnegative Borel measure on [0, 00). We shall need a multidimensional 
version of the Bernstein-Hausdorff-Widder theorem, valid for arbitrary cones. Such 
a result was proven by Choquet 

Theorem 2.2 (Bernstein— Hausdorff—Widder— Choquet theorem) Let V be a 

finite- dimensional real vector space, let C be an open convex cone in V, and let C* = 
{£ G V*: (£, x) > for all x G C} be the closed dual cone. Then a function f: C — > R 
is completely monotone if and only if there exists a positive measure /1 on C* satisfying 

f(x) = J e-fc*> dn(£) . (2.3) 
c* 

In this case, n{C*) = f{0 + ); in particular, [l is finite if and only if f is bounded. 

In particular, if f is completely monotone on C , then it is extendible [using \2. 3\) J 
to an analytic function on the complex tube C + iV satisfying 

\D Ul ...D u J(x + iy)\ < (-l) k D ul ...D u J(x) (2.4) 

for allk>Q,xEC,y^V and u 1; . . . , G C . 

15 See also Nussbaum [108] . Devinatz and Nussbaum g7], Hirsch [HI Section VII. 2], Glockner 
[6i Section 16] and Thomas [137] . 
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Remarks. 1. Since C is nonempty and open, it is not hard to see that {0} 
implies (£, x) > for all x G C. It then follows from (I2.3P that either 

(a) n is supported on {0}, in which case / is constant, 
or else 

(b) we have the strict inequality 

(-l) fc £> Ul ...£> Ufc /(z) > (2.5) 

for all k > 0, all u 1; . . . , G C, and all x G C. 

Furthermore, for Ui, . . . , in the closure of C, the left-hand side of (12. 5ft is either 
strictly positive for all x G C or else identically zero on C. Of course, the latter case 
can occur: e.g. f(xi,x 2 ) = e~ Xl on (0, oo) 2 . 

2. In our definition of complete monotonicity, the cone C plays two distinct roles: 
it is the domain on which / is defined, and it provides the direction vectors Uj for 
which the inequalities (12. ip hold. Choquet [12] elegantly separates these roles, and 
considers functions on an arbitrary open set Q C V that are completely monotone 
with respect to the cone C. He then proves the integral representation (12. 3p under 
the hypothesis fl + CCfi. This is a beautiful generalization, but we shall not need 
it. □ 

By virtue of Theorem 12.21 one way to test a function / for complete monotonicity 
is to compute its inverse Laplace transform and ask whether it is nonnegative and 
supported on C*. Of course, this procedure is not necessarily well-defined, because the 
inverse Laplace transform need not exist; moreover, if it does exist, it may need to be 
understood as a distribution in the sense of Schwartz [126J rather than as a pointwise- 
defined function. But we can say this: If /: C — > M is the Laplace transform of a 
distribution T onV*, then / is completely monotone if and only if T is positive (hence 
a positive measure) and supported on C*. This follows from the uniqueness of inverse 
Laplace transforms on the space D'(IR n ) of distributions |126[ p. 306]. Note that the 
complete monotonicity of / can fail either because T fails to be positive or because 
T fails to be supported on C*. In the former case, we can conclude that / is not 
completely monotone on any nonempty open convex subcone C C C. In the latter 
case, T might possibly be supported on some larger proper cone; but if it isn't (e.g. 
if the smallest convex cone containing the support of T is all of V*), then once again 
we can conclude that / is not completely monotone on any nonempty open convex 
subcone C C C . And finally, if / is not the Laplace transform of any distribution on 
V*, then it is certainly not the Laplace transform of a positive measure, and hence is 
not completely monotone on any nonempty open convex subcone C C C. 

In the applications to be made in this paper, the functions / will typically be of 
the form F~^, where F is a function (usually a polynomial) that is strictly positive 
on the cone C and has an analytic continuation to the tube C + iV (for polynomials 
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this latter condition of course holds trivially) . The following corollary of Theorem 12.21 
shows that the complete monotonicity of F~° on the real cone C implies the absence 
of zeros of F in the complex tube C + iV: 

Corollary 2.3 Let V be a finite- dimensional real vector space and let C be an open 
convex cone in V. Let F be an analytic function on the tube C + iV that is real and 
strictly positive on C. If F^ 13 is completely monotone on C for at least one (3 > 0, 
then F is nonvanishing on C + iV . [In particular, when V = W n and C = (0, oo) n ; 
the function F has the half-plane property.] 

Proof. Suppose that G = F~^ is completely monotone on C; then by Theorem 12.21 
it has an analytic continuation to C + iV (call it also G). Now suppose that S = {z G 
C + iV: F(z) = 0} is nonempty. Choose a simply connected domain D C (C + iV)\S 
such that D(lC ^ and 5n5/0@ Then H = F~P is a well-defined analytic 
function on D (we take the branch that is real and positive on D n C). On the other 
hand, H coincides with G on the real environment D PI C, so it must coincide with 
G everywhere in D. But lim = +oo for all z e D fl S, which contradicts the 

analyticity of G on C + iV . □ 

Remarks. 1. It also follows that the analytic function G = F~^ defined on C + iV 
is nonvanishing there. 

2. The hypothesis that F have an analytic continuation to C + iV is essential; it 
cannot be derived as a consequence of the complete monotonicity of F~^ on C. To 
see this, take V — R and C = (0, oo) and consider F(x) = (1 + |e _:c ) _1 ^ with any 
0>O. 

3. The converse of Corollary 12.31 is easily seen to be false: for instance, the 
univariate polynomial P(x) = 1 + x 2 has the half-plane property (i.e. is nonvanishing 
for Rex > 0), but P~ 13 is not completely monotone on (0, oo) for any > 0. The 
same holds for the bivariate multiaffine polynomial P(xi,Xz) = 1 + x±X2- So the 
complete monotonicity of P _/3 for some /3 > is strictly stronger than the half-plane 
property. It would be interesting to know whether similar counterexamples can be 
found if P is required to be homogeneous, or homogeneous and multiaffine. □ 

In this paper we will typically consider a polynomial P that is strictly positive 
on an open convex cone C, and we will ask for which values of ft the function P _/3 
is completely monotone. We begin with a trivial observation: If P is a nonconstant 

16 For instance, let Q, be a simply connected open subset of C whose closure is a compact subset 
of C and which satisfies (O + iV) nS/0;fka norm || • || on V; and let 

R = inf ||y|| . 

x e n, y g v 

x + iy 6 S 

By compactness we must have R > (for otherwise we would have £1 D S ^ 0, contrary to the 
hypothesis that F > on C). Now take D = Q + iBr, where B R = {y e V: \\y\\ < R}. 
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polynomial, then P~" cannot be completely monotone on any nonempty open convex 
cone for any < (because P grows at infinity in all directions except at most a 
variety of codimension 1); and P _/3 is trivially completely monotone for (3 = 0. So 
we can restrict attention to /3 > 0. 

Given a function F: C — >■ (0, oo) — for instance, a polynomial — we can ask 
about the set 

&f = {P > 0: is completely monotone onC} . (2.6) 

Clearly Bp is a closed additive subset of (0, oo). In particular, we either have Bp C 
[e, oo) for some e > or else Bp = (0,oo). The following easy lemma [TBI ED] 
characterizes the latter case: 

Lemma 2.4 Let V be a finite- dimensional real vector space, let C be an open convex 
cone in V, and let F: C — > (0, oo). Then the following are equivalent: 

(a) F~^ is completely monotone on C for all (3 > 0. 

(b) F~^ i is completely monotone on C for a sequence {Pi} of strictly positive num- 
bers converging to 0. 

(c) — logF is conditionally completely monotone on C. 

Proof, (a) => (b) is trivial, and (b) => (c) follows from 

F-P-l 

-logF = lim (2.7) 

P\a (3 

and its derivatives with respect to x. Finally, (c) => (a) follows from F~^ = 
exp(— (3 logF) and the fact that exp is absolutely monotone (i.e. has all derivatives 
nonnegative) on (— oo, oo). □ 

Already for C = (0, oo) it seems to be a difficult problem to characterize in a 
useful way the functions F described in Lemma [2.41 or even the subclass consisting 
of polynomials For polynomials P(x) = ]^[(1 + x/xj), a necessary condition 

from Corollary 12.31 is that P have the half-plane property, i.e. ReXj > for all i. A 
sufficient condition is that all Xi be real and positive; and for quadratic polynomials 
this condition is necessary as well. But already for quartic polynomials the situation 
becomes more complicated: for instance, we can take x\ = a + bi, X2 = a — hi, 
x 3 = x 4 = c with < c < a and 6 6 I, and it is not hard to see that — log P is 
conditionally completely monotone on (0, oo)@ 

It also seems to be a difficult problem to characterize the closed additive subsets 
S C (0, oo) that can arise as S = Bp. 

17 Functions / = F^ 1 for which f 13 is completely monotone for all (3 > are sometimes called 
logarithmically completely monotone 20 . 

18 The function — log P is conditionally completely monotone on (0, oo) if and only if (log P)' = 
J~l(x + Xi)^ 1 is completely monotone on (0,oo); and this happens if and only if its inverse Laplace 

i 

transform, which is g(t) = ^2 e~ tXi , is nonnegative on [0, oo). 
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Example 2.5 Fix a > 0, and consider F(x) = (1 + ae' x )~ 1 . Then the function 
Fix)' 13 = (1 + ae~ x Y is obviously completely monotone on (0, oo) whenever G 
{0, 1, 2, 3, . . .}. On the other hand, if /3 ^ {0, 1, 2, 3, . . .} we claim that F _/3 is not 
completely monotone. Indeed, for < a < 1 the convergent binomial expansion 




shows that F~P is the Laplace transform of the signed measure a fe { l D$k, which 

is nonnegative if and only if (3 G {0, 1, 2, 3, . . .}. On the other hand, for a > 1 the 
function F~@ has singularities in the right half-plane at x — log a ± in whenever 
f3 ^ {0, 1, 2, 3, . . .}, so it is not the Laplace transform of any distribution. □ 



Example 2.6 It is an interesting problem [T(H [T ^ 1601 |6"H [T0T| I149j to determine the 
pairs (/i, A) G M 2 for which the function 

F,,x(x) = x~»(x 2 + l)- x (2.9) 

is completely monotone on (0, oo). It is easy to show that there is a function /x*(A) 
such that F^\ is completely monotone if and only if fi > /x*(A); furthermore, the 
function /i* is subadditive. The state of current knowledge about /i* seems to be: 

/i*(A) = -2A for A < (2.10a) 

A < //*(A) < min(2A, 1) for < A < 1 (2.10b) 

/i*(A) = 2A + o(A) for A 1 (2.10c) 

^(A) = A for A > 1 (2.10d) 

It seems to be an open problem even to prove that //* is continuous. □ 



3 Constructions: Series and parallel connection 

In this section we discuss some general constructions by which new polynomials 
P with P^ 13 completely monotone can be obtained from old ones. In the situations 
we have in mind, the vector space V decomposes as a direct sum V = V\ © V2 and the 
cone C is a product cone C = C\ x C2 (with C\ C V\ and Ci C V2). Since we shall 
be using the letters A, B, C, D in this section to denote functions, we shall write our 
cones as C. 

Let us begin with a trivial fact: a function f(x,y) that is completely monotone 
on C\ x C2 can be specialized by fixing y to a specific point in C2, and the resulting 
function will be completely monotone on C\. In particular, this fixed value can then 
be taken to zero or infinity, and if the limit exists — possibly with some rescaling - 
then the limiting function is also completely monotone on C±. Rather than stating a 
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general theorem of this kind, let us just give the special case that we will need, which 
concerns functions of the form 

f(x,y) = [A(x) + B(x)y}~P (3.1) 

with V 2 = R and C 2 = (0, oo). 

Lemma 3.1 Let V be a finite- dimensional real vector space and let C be an open 
convex cone in V . Fix (3 > 0, and let A,B: C — > (0, oo). If (A + By)' 13 is completely 
monotone on C x (0, oo), then A' 13 and B~@ are completely monotone on C. 



Proof. Restrict to fixed y G (0, oo) and then take y I 0; this proves that A' 13 is 
completely monotone. Restrict to fixed y G (0, oo), multiply by y 13 and then take 
y t oo; this proves that B' 13 is completely monotone. □ 

As we shall see later, this trivial lemma is an analytic version of deletion {y — > 0) or 
contraction (y — > oo) for graphs or matroids. 

Let us also observe a simple but important fact about complete monotonicity for 
functions defined on a product cone C\ x C 2 : 

Lemma 3.2 For i = 1,2, let Vi be a finite- dimensional real vector space and let Ci 
be an open convex cone in Vi. Let f: C\ xC 2 —> K. Then the following are equivalent: 

(a) f is completely monotone on C\ x C%. 

(b) For all k > 0, ally G C 2 , and all choices of vectors Ui, . . . , Ufc G C 2 , the function 

Qk 

Fk,y,u 1 ,...,u k (x) = (-l) fc ^T f(x, y + hut + . . . + t k u k ) 

ot 1 ...ot k tl =...=t k =o 

(3.2) 

is completely monotone on C\. 

(c) For all k > 0, all y G C 2 , and all choices of vectors Ui, . . . ,u k e C 2 , there exists 
a positive measure fi k ,y )U1 ,...,u k on C{ such that 

F k , y , Ul ,..., Uk (x) = J< " r///,,,,,, ,,(/). (3.3) 



In particular, when Vi — R and C 2 = (0, oo), (b) reduces to the statement that the 
functions F kjy (x) = (—l) k d k f/dy k are completely monotone on C\ for all k > and 
all y > 0, and (c) reduces analogously. 

Proof of Lemma I3~2l (a) ^=^> (b) is a trivial consequence of the definition (12. lft , 
while (b) <^=^ (c) follows immediately from the the Bernstein-Hausdorff-Widder- 
Choquet theorem (Theorem 12. 2p . □ 
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Statement (c) can be rephrased loosely as saying that the inverse Laplace trans- 
form of f(x,y) with respect to x is a completely monotone function of y e C 2 . (To 
make this more precise, one should add the same qualifications as in the paragraph 
after Theorem O) 

One important application of this lemma concerns functions of the form (13. ip : 

Lemma 3.3 Let V be a finite-dimensional real vector space and let C be an open 
convex cone in V . Fix (3 > 0, and let A,B: C — > (0, oo). Then the following are 
equivalent: 

(a) (A + By)' 13 is completely monotone on C x (0, oo). 

(b) B~P exp(—tA/B) is completely monotone on C for all t > 0. 

(c) B K (A + zB)~(P +K ) is completely monotone on C for all k > —(3 and all z > 0. 

(d) B k (A + zB)~^ +k ' is completely monotone on C for all integers k > and all 
z>0. 



Proof. We have the Laplace-transform formula 

oo 

(A + By)'? = J dt e~ ty B^ exp(-tA/B) . (3.4) 
o 

Therefore, Lemma [3.2( a) <^=^ (c) with C\ = (0, oo) and C2 = C proves the equiva- 
lence of (a) and (b). 

Now assume that B^ 13 exp(— tA/B) is completely monotone on C for all t > 0. 
Then we can multiply by e~ 2 *t p_1 /r(p) for any p > and integrate over t G (0, 00), 
and the result will be completely monotone. This (together with a trivial evaluation 
at t = to handle p — 0) shows that (b) =^ (c). 

(c) =>■ (d) is trivial. 

The equivalence of (a) and (d) follows from Lemma l3T2T a) <^=^ (b), used with 
d = C and C 2 = (0,oo). □ 



Corollary 3.4 LetV be a finite-dimensional real vector space, letC be an open convex 
cone in V , and let A, B: C — > (0, 00). Define 

= {(3 > 0: B^ 13 is completely monotone on C} (3.5) 
&A+B y = {P > 0: (A + By)' 13 is completely monotone on C x (0, 00)} (3.6) 

Then B A+By + B B C B A+By . 

In particular, if B B = (0, 00), then Ba+b v is either the empty set or all of (0, 00) 
or a closed interval [(3 , 00) with O > 0. 
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Proof. This follows immediately from Lemma [3.3( a) •<=>■ (b): for if ft G Ba+bv and 
A G Bb, then exp(— tA/B) and B~ x are both completely monotone on C, hence 
so is their product, hence ft + A G Ba+b v - D 



Lemma [3.31 leads to the following extremely important result, which (as we shall 
see later) is an analytic version of parallel connection for graphs or matroids: 

Proposition 3.5 Let V be a finite-dimensional real vector space and let C be an 

open convex cone in V. Fix (3 > 0, and let A, B,C, D: C — > (0, oo). Suppose that 
(A + By)' 13 and (C + Dy)' 13 are completely monotone on C x (0, oo). Then the same 
is true of (AD + BC + BDy)~@. 

Proof. By Lemma l3.3[ B~@ exp(— tA/B) and D~@ exp(— tC/D) are completely mono- 
tone on C for all t > 0. Hence the same is true of their product, which is 
(BD)' 13 exp[— t(AD + BC) / (BD)]. But then using Lemma 1331 again (this time in the 
reverse direction), we conclude that (AD + BC + BDy)' 13 is completely monotone 
onCx(0,oo). □ 

We also have an analytic version of series connection for graphs or matroids, but 
only for ft > 1/2: 

Proposition 3.6 Let V be a finite-dimensional real vector space and let C be an 

open convex cone in V. Fix ft > 1/2, and let A, B,C, D: C — > (0, oo). Suppose that 
(A + By)' 13 and (C + Dy)' 13 are completely monotone onC x (0, oo). Then the same 
is true of [AC + (AD + BC)y]~ p . 



To prove Proposition 13.61 we begin with a lemma that we think is of independent 
interest; both the sufficiency and the necessity will play important roles for us. 



Lemma 3.7 For ft G R and A > 0, the function 



(u + v) 13 exp^— A 



uv 



u + v 



(3.7) 



is completely monotone on (0, oo) 2 if and only if ft > 1/2. 

In particular, for ft > 1/2 there exists a positive measure on [0, oo) 2 such that 



(u + v)-? 



exp 



-A- 



uv 



u + v 



-t\U—tlV 



(3.8) 



[0,oo) 2 



PROOF. "If": Since (u + v)~^~^ is completely monotone when ft > 1/2, it suffices 
to prove the complete monotonicity for ft = 1/2. But this follows immediately from 
the identity 



(u+v) 1//2 exp 



-A- 



uv 
u + v 



—= I ds exp 



S + -) u 



— i 
2 / 



(3.9) 
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which is easily verified by completing the square in the Gaussian integral. The state- 
ment about the measure /i^a then follows from the Bernstein-Hausdorff-Widder- 
Choquet theorem (Theorem 12. 2p . 

"Only if": If Fp t \ is completely monotone, then so is Fp> t x for all f3' > (3; so it 
suffices to prove the failure of complete monotonicity for < /3 < 1/2. Now, by 
Lemma [3.21 Fp t \ is completely monotone on (0, oo) 2 if and only if the functions 

gk 



(3.10) 



are completely monotone on (0, oo) for all k > and all v > 0, or equivalently if their 
inverse Laplace transforms with respect to u, 



Qk 



(3.11) 



(see [5*2] p. 245, eq. 5.6(35)]), are nonnegative for all k > and all t, v > 0. For k = 
this manifestly holds for all (3 > 0; but let us now show that for k — 1 it holds only 
for p > 1/2. We have 



t + X + 







Ip-i(2vVXt) - 2VXtI' /3 _ l (2vVXt) 

(3.12) 



and we need the term in square brackets to be nonnegative for all t, v > 0. Write 
x = 2v\[Xt and eliminate t> in favor of x; we need 



t + 2VAt 



/3 



X 



+ A > 



I/3-l(x 

for all t, x > 0. This quadratic in \/t is nonnegative for all t > if and only if 

0-1 



(3.13) 



But using the large- a; asymptotic expansion 



> -1 . 



(3.14) 



J-log/^Or) = 1 - i- + 0(l/x 2 ) 



we see that 



/3-1 4_i (^) 



-1 + 



/3 



+ 0(l/a; 2 



X Ip_i(x) X 

which is < — 1 for all sufficiently large x whenever (3 < 1/2. □ 



(3.15) 



(3.16) 



Remarks. 1. It is obvious by rescaling of u and v that, for any given {3, the functions 
Fp^x are either completely monotone for all A > or for none. 
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2. The appeal to the Bernstein-Hausdorff-Widder-Choquet theorem can be 
avoided: for (3 = 1/2, the integral representation (13. 9 p already provides the desired 
measure Hi/2,\] and for (3 > 1/2, (13.91) together with 

{u + vY^h = Jdt-L^e-^) (3.17) 
o 2 

represents /ip^ as the convolution of two positive measures. Indeed, multiplying (13. 9p 
by (I3.17p . one obtains after a straightforward change of variables the explicit formula 

(AX) 1 - 13 a 3 

todtuh) = T ^ )T ^_i ) P(ti,t 2 ,xf-2 X (tut 2 ,x) (3.18) 



where 



and 



P(£i,£ 2 ,A) = 2(tit 2 + Xti + Xt 2 ) - (t{ + t 2 + A ) (3.19) 



,1 if t u t 2 > and — Vh\ < VX < + ^T 2 , . 

X{ti,t 2 ,A) — < [6.Z0) 



otherwise 

The constraint x(ti,h,X) ^ states simply that y/t[, y/t^, \/X form the sides of a 
triangle; and P(ti,t 2 ,X) is precisely_16 times the square of the area of this triangle 
(Heron's formula [HI Section 3.2] )o In view of these explicit formulae, the proof of 
Lemma 13.71 is in fact completely elementary. 

3. It would be interesting to know whether Lemma 13.71 can be generalized to 
other ratios of elementary symmetric polynomials, e.g. E~P exp(— XE r+ i jn /E rjn ). By 
Lemma I3T31 this would determine the complete monotonicity of E~+ ln+l . □ 



Proof of Proposition GLU By LemmaES! B^ 13 exp(—tiA/B) and D~P exp(-t 2 C/D) 
are completely monotone on C for all t\,t 2 > 0. Hence the same is true of their 
product for any choice of £1,^2 > 0. We now use the identity (I3.8p . multiplied on 
both sides by (BD)" 13 , with u = A/B and v = C/D. This shows that (AD + 
BCy l3 exp[-XAC/(AD + BC)] is completely monotone on C for all A > 0. But 
then using Lemma 13.31 again (this time in the reverse direction), we conclude that 
[AC + (AD + BC)y]~P is completely monotone on C x (0, 00). □ 

Using Lemma 13.71 we can also show that the spanning-tree polynomial of the 3- 
cycle — or equivalently, the elementary symmetric polynomial of degree 2 in three 
variables — has the property that P _/3 is completely monotone on (0, oo) 3 if and only 
ifP> 1/2: 

19 It is curious that similar expressions, involving the area of a triangle in terms of its sides, arise 
also in Sonine's integral for the product of three Bessel functions [1441 p. 411, eq. 13.46(3)] ITJ p. 36, 
eq. (4.39) and p. 40] — a formula that Szego [135] employed in one version of his nonnegativity proof 
for (|1.1[) in the case n = 3. Probably this is not an accident; it would be interesting to understand 
the precise connection between Sonine's formula and (|3.18|) [cf. also (|3.21[) below]. 
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Proposition 3.8 The function F(xi,X2,Xs) = (x%X2 + X1X3 + X2X3) 13 is completely 
monotone on (0, oo) 3 if and only if (3 = or (3 > 1/2. 



Proof. As always, F is completely monotone for (3 = and not completely monotone 
for (3 < 0, so it suffices to consider (3 > 0. Using Lemma 13.31 with A = X1X2, 
B = x\ + X2, y = £3, we see that F is completely monotone on (0, oo) 3 if and only 
if the function F@ t x defined in (13.71) is completely monotone for all A > 0. But by 
Lemma [3 .7\ this occurs if and only if > 1/2. □ 



Remarks. 1. We shall later give two further independent proofs of Proposition 13.81 
one based on harmonic analysis on the cone of positive-definite mxm real symmetric 
matrices specialized to m = 2 [Corollaries II .51 and 16.9( a)] . and one based on harmonic 
analysis on the Lorentz cone in R n specialized to n = 3 [Corollary 11.10} proved in 
Section [7J. The point here is that the Jordan algebra Sym(2, R) ~ 1 x R 2 can be 
viewed as a member of two different families of Jordan algebras: Sym(m, R) and 
R x W 1 - 1 [56, p. 98]. 

2. Proposition 13.81 implies that the property stated in Proposition 13.61 does not 
hold for < (3 < 1/2. Indeed, it suffices to take C = (0, oo) 2 and A = xi, B = 1, 
C = X2, D = 1, leading to the function (x±X2 + x%y + X2y)~^ ■ 

3. Proposition 13.81 is, of course, just Theorem 11.11 restricted to the 3-cycle G = 
K 3 . In particular it implies Szego's |135] result (except the strict positivity) for the 
polynomial (II .ip in the special case n = 3. 

4. Combining (J33D with (ETSjl /f lBTBj) . we obtain for > 1/2 the formula 

+ = r(i)r ( g-i)r(g) x 

00 00 00 

y dtx J dt 2 j dfc 3e -*i*i-^-*a*a Pfat^tffi , (3.21) 


which provides an explicit elementary proof of the direct ( "if" ) half of Proposition ^. 81 
See also Remark 1 in Section 16731 for an alternate derivation of (I3.2ip . and see Corol- 
lary [722 for a generalization from to E 2 , n - CJ 

The following generalization of Lemma I3~71 is an analytic version of series extension 
of a single edge: 

Lemma 3.9 Let V be a finite- dimensional real vector space and let C be an open 
convex cone in V. Let f be completely monotone on C x (0, 00), and let (3 > 1/2. 
Then the function 

f*(x,y,z) = ( y + z )-f>f( x ,J!L-) (3.22) 
is completely monotone on C x (0, oo) 2 . 
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Proof. By the Bernstein-Hausdorff-Widder-Choquet theorem (Theorem 12 .2 j) and 
linearity, it suffices to prove the lemma for f(x,y) = exp(—(£,x) — Xy) with £ G C* 
and A > 0. The variable x now simply goes for the ride, so that the claim follows 
immediately from Lemma 13.71 □ 



4 Application to graphs and matroids 
4.1 Graphs 

Let G = (V, E) be a finite undirected graph with vertex set V and edge set E; in 
this paper all graphs are allowed to have loops and multiple edges unless explicitly 
stated otherwise. Now let x = {x e } e& E be a family of indeterminates indexed by the 
edges of G. If G is a connected graph, we denote by Tq(x) the generating polynomial 
of spanning trees in G, namely 

r G (x) = n 

TGT(G) eGT 

where T(G) denotes the family of edge sets of spanning trees in G. If G is discon- 
nected, we define T G (x) to be the product of the spanning-tree polynomials of its 
connected components. Otherwise put, T G (x) is in all cases the generating polyno- 
mial of maximal spanning forests in G. This is a slightly nonstandard definition (the 
usual definition would put T G = if G is disconnected), but it is convenient for our 
purposes and is natural from a matroidal point of view (see below). In order to avoid 
any possible misunderstanding, we have inserted in Theorems 11.11 and 11.11 and Corol- 
lary [L8] the word "connected", so that the claims made in the Introduction will be 
true on either interpretation of Tq- Please note that, in our definition, T G is always 
strictly positive on (0, oo) £ , because the set of maximal spanning forests is nonempty. 
Note also that, on either definition of T G , loops in G (if any) play no role in T G . And 
it goes without saying that T G is a multiaffine polynomial, i.e. of degree at most 1 in 
each x e separately. 

If e is an edge of G, the spanning-tree polynomial of G can be related to that of 
the deletion G\e and the contraction G/e: 

(T G \ e (x^ e ) + x e T G / e (x^ e ) if e is neither a bridge nor a loop 
x e T G \ e (x^ e ) = x e T G/e (x^ e ) if e is a bridge 
7 G \e(x^ e ) = T G/e (x^ e ) if e is a loop 

(4.2) 

where x^ e denotes {^/}/ee\{ e }- The fact that T G \ e equals T G / e (rather than equalling 
zero) when e is a bridge is a consequence of our peculiar definition of T G . 

Now let us take an analytic point of view, so that the indeterminates x e will be 
interpreted as real or complex variables. 

Definition 4.1 For each (3 > 0, we denote by Qp the class of graphs G for which 
Tq^ is completely monotone on (0, oo)" B . 
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The naturality of the classes Qp is illustrated by the following easy but fundamental 
result: 

Proposition 4.2 Each class Qp is closed under taking minors (i.e. under deletion 
and contraction of edges and deletion of isolated vertices ), under disjoint unions, and 
under gluing at a cut vertex. 

Proof. Deletion of a non-bridge edge e corresponds to taking x e 4- 0. Contraction of 
a non-loop edge e corresponds to dividing by x e and taking x e ~\ +oo. Both of these 
operations preserve complete monotonicity. Deletion of a bridge has the same effect 
as contracting it, in our peculiar definition of T G . Contraction of a loop is equivalent 
to deleting it (but loops play no role in T G anyway). Isolated vertices play no role in 
Tq. This proves closure under taking minors. 

If G is obtained from G\ and G 2 either by disjoint union or by gluing at a cut 
vertex, then T G = T Gl T G2 (on disjoint sets of variables) in our definition of T G ; this 
again preserves complete monotonicity. □ 

Proposition 14.21 illustrates the principal reason for allowing arbitrary constants 
c > (rather than just c = 1) in Theorem 11.11 and subsequent results: it leads 
to a minor-closed class of graphs. This, in turn, allows for characterizations that 
are necessary as well as sufficient. A similar situation arises in studying the negative- 
correlation property for a randomly chosen basis of a matroid. If only the "uniformly- 
at-random" situation is considered (i.e., element weights x = 1), then the resulting 
class of matroids is not minor-closed, and closure under minors has to be added by 
hand, leading to the class of so-called balanced matroids [59J. But it then turns out 
that the class of balanced matroids is not closed under taking 2-sums [U]. If, by 
contrast, one demands negative correlation for all choices of element weights x > 
0, then the resulting class — the so-called Rayleigh matroids - - is automatically 
closed under taking minors (by the same x e — > and x e — > oo argument as in 
Proposition I4.2p . Moreover, it turns out to be closed under 2-sums as well [4"Tj . 

The very important property of closure under 2-sums holds also in our context. 
To see this, let us first recall the definitions of parallel connection, series connection 
and 2-sum of graphs |110l Section 7.1], and work out how T G transforms under these 
operations. 

For i — 1, 2, let Gj = (V*, Ej) be a graph and let be an edge of Gf, it is convenient 
(though not absolutely necessary) to assume that e, is neither a loop nor a bridge 
in G{. Let us furthermore choose an orientation e*j = x^y\ for the edge e%. (To avoid 
notational ambiguity, it helps to assume that the sets V\, V 2 , Ei, E 2 are all disjoint.) 
Then the parallel connection of (Gi,e*i) with (G 2 ,e 2 ) is the graph (Gx, ei)||(Cr 2 , e 2 ) 
obtained from the disjoint union Gi UG 2 by identifying x\ with x 2 , y\ with y 2 , and e± 
with e 2 . [Equivalently, it is obtained from the disjoint union (Gi \ e\) U {G 2 \ e 2 ) by 
identifying x\ with x 2 (call the new vertex x), y\ with y 2 (call the new vertex y), and 
then adding a new edge e from x to y.} The series connection of (Gi, e*i) with (G 2 , e 2 ) 
is the graph e\) ix (G 2 , e 2 ) obtained from the disjoint union (Gi\ei)U(G 2 \e 2 ) by 
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identifying x\ with x 2 and adding a new edge e from y\ to y 2 . The 2-sum of (Ci, ei) 
with (G2,e 2 ) is the graph (G^ei) © 2 iG 2l e 2 ) obtained from the parallel connection 
(Gi, e*i) 11(^2, e 2 ) by deleting the edge e that arose from identifying e x with e 2 , or 
equivalently from the series connection e*i) M (G 2 , e 2 ) by contracting the edge e. 

To calculate the spanning-tree polynomial of a parallel connection, series connec- 
tion or 2-sum, it is convenient to change slightly the notation and suppose that e\ and 
e 2 have already been identified (let us call this common edge e), so that E 1 nE 2 = {e}. 
It is then not difficult to see [1101 Proposition 7.1.13] that the spanning-tree polyno- 
mial of a parallel connection Gi|| e (j2 is given by 

^Gi|| e G 2 ( x ) = T Gl \ e T G2 / e + T Gl/e T G2 \ e + x e T Gl / e T G2 / e , (4.3) 

while that of a series connection G± cxi e G 2 is 

T Gltxi£G2 {x) = T Gl \ e T G2 \ e + x e T Gl \ e T G2 / e + x e T Gl / e T G2 \ e . (4.4) 

(All the spanning-tree polynomials on the right-hand sides are of course evaluated at 
x^ e .) The spanning-tree polynomial of a 2-sum G\ ©2,e G 2 is therefore 

^Gi®2, e G 2 (x) = T Gl \ e T G2 / e + T Gl/e T G2 \ e . (4.5) 

Proposition 4.3 Each class Qp is closed under parallel connection and under 2- 
sums. 

Proof. Closure under parallel connection is an immediate consequence of Proposi- 
tion 13.51 and the formula H4.3[) for parallel connection. Since the 2-sum is obtained 
from the parallel connection by deletion, closure under 2-sum then follows from Propo- 
sition S2J □ 

Proposition 4.4 The class Qp is closed under series connection for f3 > 1/2 but not 
forO< < 1/2. 

Proof. For > 1/2, closure under series connection is an immediate consequence 
of Proposition 13.61 and the formula (14.41) for series connection. For < < 1/2, 
non-closure under series connection follows immediately from the observation that 
the 2-cycle C 2 = K { 2 2) (a pair of vertices connected by two parallel edges) lies in Qp 
for all > 0, but the series connection of a 2-cycle with another 2-cycle is a 3-cycle, 
which lies in Qp only for (3 > | (by Proposition 13.81) . □ 

Remarks. 1. Unlike the situation for the half-plane and Rayleigh properties, the 
classes Qp are not in general closed under duality for planar graphs. For instance, the 
graph C* = Kf ] (a pair of vertices connected by three parallel edges) lies in Qp for 
all (3 > 0; but its dual C 3 lies in Qp only for (3 > ~ (by Proposition 13. 8j) . 

However, the class Qp is duality-invariant for (3 G (|, 1), as it consists of all series- 
parallel graphs (Theorem 15.31) . And since Qp = all graphs for (3 G {|, 1, |, . . .}, these 
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classes QpHV (where V = planar graphs) are also duality- invariant. We do not know 
whether Gp D V is duality-invariant for (3 G (1, oo) \ {|, 2, |, . . .}. In any case, the 
duality question is most naturally posed for matroids rather than for graphs. 

2. Since Gp is closed under 0-sums (disjoint unions), 1-sums (gluing at a cut 
vertex) and 2-sums (essentially gluing at an edge), it is natural to ask whether it is 
also closed under 3-sums (gluing along triangles). We do not know the answer. In 
particular, K4 G Gp for all > 1 by Corollary ll.5[ but as noted in the discussion 
after Problem 11.131 . we do not know whether K 5 — e = K4 © 3 K4 belongs to Gp for 

0e (!,§)■ □ 

It is an well-known (and easy) result that any minor-closed class of graphs is of 
the form 

Ex (J 7 ) = {G: G does not contain any minor from J 7 } (4-6) 

for some family T of "excluded minors"; indeed, the minimal choice of T consists 
of those graphs that do not belong to the class in question but whose proper minors 
all do belong to the class. (Here we consider isomorphic graphs to be identical, or 
alternatively take only one representative from each isomorphism class.) 

In one of the deepest and most difficult theorems of graph theory, Roberston and 
Seymour |122] sharpened this result by proving that any minor-closed class of graphs 
is of the form Ex(J r ) for some finite family J 7 . Therefore, each of our classes Gp can 
be characterized by a finite family of excluded minors. One of the goals of this paper 
- alas, incompletely achieved — is to determine these excluded minors. 

4.2 Matroids 

The foregoing considerations have an immediate generalization to matroids. (Read- 
ers unfamiliar with matroids can skip this section without loss of logical continuity.) 
Let M be a matroid with ground set E, and let x = {x e } e( zE be a family of indeter- 
minates indexed by the elements of M. We denote by Bm(x) the basis generating 
polynomial of M, namely 

%(x) = n ii 1 ' ( 4 - r ) 

BeB(M) e€B 

where B(M) denotes the family of bases of M. Please note that loops in M (if any) 
play no role in Bm- Note also that if M is the graphic matroid M(G) associated to a 
graph G, we have 5m(g)(x) = T G (x.). This identity would not hold for disconnected 
G if we had taken the standard definition of Tq. 

Definition 4.5 For each (3 > 0, we denote by Aip the class of matroids M for which 
B^f is completely monotone on (0,oo) E . 

Once again we have: 

Proposition 4.6 Each class M.p is closed under taking minors (i.e. deletion and 
contraction) and direct sums. 
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Proof. The proof is identical to that of Proposition 14.21 if one substitutes "element" 
for "edge" , "coloop" for "bridge" , and "direct sum" for either form of union. □ 

We refer to |110[ Section 7.1] for the definitions of parallel connection, series 
connection and 2-sum of matroids, which generalize those for graphs. The upshot [1101 
Proposition 7.1.13] is that the formulae (I4.3p — (I4.5p for the spanning-tree polynomials 
of graphs extend unchanged to the basis generating polynomials of matroids. We 
therefore have: 

Proposition 4.7 Each class Aip is closed under parallel connection and under 2- 
sums. 

Proposition 4.8 The class M.p is closed under series connection for (3 > 1/2 but 
not for < p < 1/2. 

Since each Aip is a minor-closed class, we can once again seek a characterization 
of M.p by excluded minors. However, in this case no analogue of the Robertson- 
Seymour theorem exists, so we have no a priori guarantee of finiteness of the set 
of excluded minors. Indeed, there exist minor-closed classes of matroids having an 
infinite family of excluded minors [ 1101 Exercise 6.5.5(g)]; and in fact, for any infinite 
field F, the class of F-representable matroids has infinitely many excluded minors 
[TTOl Theorem 6.5.17]. 

We suspect that the classes M.$ are not closed under duality in general. For 
instance, Corollary 11.101 shows that f/2,5 G M.p if and only if (3 > 3/2; but we suspect 
(Conjecture II. lip that G M.p if and only if (3 > 1. On the other hand, we shall 
show in the next section (Theorem 15 .4p that Aip for 1/2 < j3 < 1 consists precisely of 
the graphic matroids of series-parallel graphs — a class that is closed under duality. 

5 Series-parallel graphs (and matroids): Proof of 
Theorem II. If 

Before proving Theorem ll.il . let us prove a similar but simpler theorem concerning 
the interval < (3 < |. 

Theorem 5.1 Let G be a graph, and let (3 G (0, |). Then the following are equivalent: 

(a) G G G p . 

(b) G can be obtained from a forest by parallel extensions of edges (i.e., replacing 
an edge by several parallel edges) and additions of loops. 

(c) G has no K 3 minor. 

Moreover, these equivalent conditions imply that G G Gp for all (3' > 0. 
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Proof. The equivalence of (b) and (c) is an easy graph-theoretic exercise. 

If G is obtained from a forest by parallel extensions of edges and additions of 
loops, then T G (x) is a product of factors of the form x ei + . . .+ x ek (where ei, . . . , ejt 
are a set of parallel edges in G), so that Tq is completely monotone on (0, oo) B for 
all > 0. Therefore (b) (a). 

Conversely, Proposition 13.81 tells us that K 3 Qp for (3 G (0, \). Since ^ is a 
minor-closed class, this proves that (a) ==>• (c). □ 

Dave Wagner has pointed out to us that Theorem 15.11 extends easily to matroids: 

Theorem 5.2 Let M be a matroid, and let f3 G (0, |). Then the following are equiv- 
alent: 

(a) M G Mp. 

(b) M is the graphic matroid M(G) for a graph G that can be obtained from a forest 
by parallel extensions of edges and additions of loops. 

(c) G has no M(K 3 ) or Uia minor. 

Moreover, these equivalent conditions imply that M G A4p> for all 0' > 0. 

Proof. Tutte has proven [1101 Theorem 10.3.1] that a matroid is graphic if and only 
if it has no minor isomorphic to U 2 a, F 7 , F 7 *, M*(K 5 ) or M*(K 3>3 ). Since M(K 3 ) is 
a minor of the last four matroids on this list, the equivalence of (b) and (c) follows 
from the graphic case. 

(b) ==>■ (a) has already been proven in Theorem 15.11 

Finally, Proposition EHJtells us that M(K 3 ) <£ Mp for /3 G (0, |), and CorollaryOl 
or 11.101 (to be proven in Sections [6] and [3, respectively) tells us that Via -M-p for 
/3 G (0, 1). Since M.$ is a minor-closed class, this proves that (a) =^ (c). □ 

Let us now prove the corresponding characterization for | < (3 < 1, which is a 
rephrasing of Theorem 11.11 and concerns series-parallel graphs. Unfortunately, there 
seems to be no completely standard definition of "series-parallel graph"; a plethora 
of slightly different definitions can be found in the literature [SHI SH 11091 IHOl I3D] . 
So let us be completely precise about our own usage: we shall call a loopless graph 
series-parallel if it can be obtained from a forest by a finite sequence of series and 
parallel extensions of edges (i.e. replacing an edge by two edges in series or two 
edges in parallel). We shall call a general graph (allowing loops) series-parallel if its 
underlying loopless graph is series-parallel] 20 ! 



20 Some authors write "obtained from a tree", "obtained from K2" or "obtained from C2" in place 
of "obtained from a forest" ; in our terminology these definitions yield, respectively, all connected 
series-parallel graphs, all connected series-parallel graphs whose blocks form a path, or all 2-connected 
series-parallel graphs. See [30[ Section 11.2] for a more extensive bibliography. 
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So we need to understand how the spanning-tree polynomial 2b (x) behaves under 
series and parallel extensions of edges. Parallel extension is easy: if G' is obtained 
from G by replacing the edge e by a pair of edges e\ and e 2 in parallel, then 

T G >(-x^ e ,x ei ,x e2 ) = T G (x^ e ,x ei +x e2 ) . (5.1) 

In other words, two parallel edges with weights x ei and x e , 2 are equivalent to a single 
edge with weight x ei +x e2 . This is because the spanning trees of G' are in correspon- 
dence with the spanning trees T of G as follows: if T does not contain e, then leave 
T as is (it is a spanning tree of G'); if T does contain e, then adjoin to T \ e one but 
not both of the edges t\ and e 2 . 

Series extension is only slightly more complicated: if G' is obtained from G by 
replacing the edge e by a pair of edges e± and e 2 in series, then 

T G /(x^ e ,x ei ,x e2 ) = (x ei + x e jT G (x^ e , Xei * e2 j, (5.2) 

In other words, two series edges with weights x ei and Ou are equivalent to a single 
edge with weight x ei x e2 /(x ei +x e2 ) together with a prefactor that clears the resulting 
denominator. This is because the spanning trees of G' are in correspondence with the 
spanning trees T of G as follows: if T does not contain e, then adjoin to T \ e one 
but not both of the edges e\ and e 2 ; if T does contain e, then adjoin to T \ e both of 
the edges e x and e 2 . Since T G (x) = T G \ e (x^ e ) + a; e T G/e (x^ e ) where T G \ e (resp. T G/e ) 
counts the spanning trees of G that do not (resp. do) contain e (the latter without 
the factor x e ), this proves (15. 2p . 

Let us remark that the parallel and series laws for T G (x) are precisely the laws for 
combining electrical conductances in parallel or series. This is no accident, because as 
Kirchhoff [90] showed a century-and-a-half ago, the theory of linear electrical circuits 
can be written in terms of spanning-tree polynomials (see e.g. [3H] for a modern 
treatment). Let us also remark that the parallel and series laws for T G (x) are limiting 
cases of the parallel and series laws for the multivariate Tutte polynomial Z G (g, v), 
obtained when q — > and v is infinitesimal: see |130[ Sections 4.4-4.7] for a detailed 
explanation. 

We are now ready to state and prove the main result of this section: 

Theorem 5.3 LetG be a graph, and let (3 G (|, 1). Then the following are equivalent: 

(a) G e Q p . 

(b) G is series-parallel. 

(c) G has no K4 minor. 

Moreover, these equivalent conditions imply that G G Qp for all (3' > \. 

Please note that Theorems 15.11 and 15.31 together imply Theorem ll.il . 
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Proof of Theorem 15.31 The equivalence of (b) and (c) is a well-known graph- 
theoretic result [IHJ Exercise 7.30 and Proposition 1.7.4] (see also [50 | 1109] ). 

Now let G be a series-parallel graph: this means that G can be obtained from a 
forest by series and parallel extensions of edges and additions of loops. As shown in 
Theorem I5.1[ if G is a forest, then T^ 13 is completely monotone for all f3 > 0. Parallel 
extension (15.11) obviously preserves complete monotonicity. By Lemma 13. 9[ series 
extension (15.21) preserves complete monotonicity whenever (3 > |. Finally, additions 
of loops do not affect T G . Therefore, every series-parallel graph G belongs to the class 
Qp for all > \. 

Conversely, Corollary 16.9( a) below (whose proof in Section [6] will not use Theo- 
rem [573]) tells us that K± ^ Qp for (3 G (|, 1). Since Qp is a minor-closed class, this 
proves (a) =>■ (c). □ 

Remarks. 1. Instead of using series and parallel extension of edges [eqns. ( 15 . 1 p / ( 15. 2 p 
and Lemma 13. 9] , we could equally well have written this proof in terms of the more 
general concept of series and parallel connection of graphs [eqns. ( 14.31) / ( 14~4| and 
Propositions 13.51 and 13. 6] . 

2. It ought to be possible to give an "elementary" proof of the fact that K4 G' Qp 
for (3 G (|,1) — and more generally of the fact that the derivatives of at c = 1 
do not all have sign (— l) k — by asymptotic calculation of coefficients a la Pemantle- 
Wilson-Baryshnikov [III ESI ESI ESI US] . □ 

Once again, Dave Wagner has pointed out to us that Theorem 15.31 extends imme- 
diately to matroids: 

Theorem 5.4 Let M be a matroid, and let (3 G (|, 1). Then the following are equiv- 
alent: 

(a) M G Mp. 

(b) M is the graphic matroid M(G) [or equivalently the cographic matroid M*(G)] 
of a series-parallel graph G. 

(c) G has no MlK^) or f/2,4 minor. 

Moreover, these equivalent conditions imply that M G M.p> for all (3' > |. 

The proof is completely analogous to that of Theorem 15.21 (see also [1101 Corol- 
lary 12.2.14] for an alternative proof of (b) <^=^> (c) in this case). 

Remark. We will prove in the next section that Qo — Q1/2 — Gi — Q3/2 — ■ ■ ■ — 
all graphs. But what is the story for matroids? Does Aii/2 contain only regular 
matroids? Does A4i contain only complex-unimodular matroids? We suspect that 
the answer to this last question is no, since we suspect that U n -2,n G A4i for all ri > 2 
(Conjecture II. lip . 
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6 Determinantal polynomials: Proof of Theorems 
IT2L O and 131 



In this section we consider polynomials defined by determinants as in (jl.4p . We 
begin with some preliminary algebraic facts about such determinantal polynomials. 
After this, we turn to the analytic results that are our principal concern. We first prove 
a simple abstract version of the half-plane property for determinantal polynomials. 
Then we turn to the main topic of this section, namely, the proof of Theorems 11.2} 
PI and PI 

6.1 Algebraic preliminaries 

First, some notation: We write [m] = {1, . . . , m}. If A = (c%)™- = i is an m x m 
matrix and I,JQ [m], we denote by Ajj the submatrix of A corresponding to the 
rows I and the columns J, all kept in their original order. We write I c to denote 
the complement of / in [m]. Then e(J, J) = (— l)^e/ l +T, jeJ j j s ^he sign of the 
permutation that takes II C into JJ C . 

We begin with a simple formula for the determinant of the sum of two matri- 
ces, which ought to be found in every textbook of matrix theory but seems to be 
surprisingly little known F*1 

Lemma 6.1 Let A,B be m x m matrices with elements in a commutative ring R. 
TheiE 

det(A + B) = e(I,J)(detA IJ )(detB I c J c) . (6.1) 

I,J C [m] 
\I\ = \J\ 

Proof. Using the definition of determinant and expanding the products, we have 
det(A + B) = Y sg 11 ! 71 ") II^W II hi '^ i ') ' ( 6 ' 2 ) 

7re<S,„ IQm] iei i'ei c 

where the outermost sum runs over all permutations tt of [m]. Define now J = 
Then we can interchange the order of summation: 

det(A + B) = 5^ s S n ( 7r ) II ai7r « II hi '^ i ') ■ ( 6 ' 3 ) 

I,JC[m] TveS m iei 

\I\ = \J\ n[I] = J 



21 This formula can be found in [97l pp. 162-163, Exercise 6] and (93j pp. 221-223]. It can also 
be found — albeit in an ugly notation that obscures what is going on — in [97l pp. 145-146 and 
163-164] 125[ pp. 31-33] |119( pp. 281-282]; and in an even more obscure notation in [3j p. 102, 
item 5]. We would be grateful to readers who could supply additional references. The proof here is 
taken from [32j Lemma A.l]. 

We remark that, by the same method, one can prove a formula analogous to (|6.ip in which all 
three occurrences of determinant are replaced by permanent and the factor e(J, J) is omitted. 

22 The determinant of an empty matrix is of course defined to be 1. This makes sense in the 
present context even if the ring R lacks an identity element: the term I = J = contributes det B 
to the sum (|6.1| . while the term I = J = [m] contributes det A. 
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Suppose now that |/| = | J\ = k, and let us write I = {z'x, . . . , ik] and J = {jx, . . . ,j k } 
where the elements are written in increasing order, and likewise I c = {i' x ,... ,i' m _ k } 
and J = . . . ,j' m _ k }- Let n' G S k and n" G S m -k be the permutations defined so 
that 

7r'(a) = < — > ir(i a ) = j p (6.4a) 
7T» = /3 <— ► tt(0=J'J (6-4b) 

It is easy to see that sgn(7r) = sgn(7r') sgn(7r")e(J, J). The formula (16. ip then follows 
by using twice again the definition of determinant. □ 

The following special case is frequently useful: 

Corollary 6.2 Let A,B be m x m matrices with elements in a commutative ring R, 
with at least one of them being a diagonal matrix. Then 

&et(A + B) = (detA//)(det B/cjc) . (6.5) 

IC[m] 

Proof. If A is diagonal, then det Ajj = whenever I ^ J, and likewise for B. □ 

Remark. We will see in Appendix[A]that Corollary 16 . 21 generalizes to quaternions, 
but Lemma [6.11 does not: see Proposition IA.5I and the remark following it. □ 

Iterating Lemma 16.11 we obtain a formula for the determinant of a sum of n 
matrices: 

(n \ n 

= "£e(I,J)l[det[(A k ) Ik j k ), (6.6) 
k=l J I,J k=l 

where the sum runs over ordered partitions I = (Ji, . . . , I n ) and J = (Jx, . . . , J n ) of 
[m] into n possibly empty blocks satisfying \I k \ = \ J k \ for all k] here e(I, J) is the sign 
of the permutation taking l x l 2 ■ ■ ■ I n into Ji J 2 • • ■ J n . 

We can now say something about determinantal polynomials of the type (II. 4p . 
Recall that if A is a (not-necessarily-square) matrix with elements in a commutative 
ring R, then the (determinantal) rank of A is defined to be the largest integer r such 
that A has a nonzero r x r minor; if no such minor exists (i.e., A = 0), we say that 
A has rank 0. 

Proposition 6.3 Let A±, . . . , A n be mx m matrices with elements in a commutative 
ring R, and let x±, . . . , x n be indeterminates. Then 

P(xi,...,x n ) = det^Y^XiA^j (6.7) 

is a homogeneous polynomial of degree m with coefficients in R. Furthermore, the 
degree of P in the variable Xi is < rank(Aj). [In particular, if each is of rank at 
most 1, then P is multiaffine.] 
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Proof. Both assertions about P are immediate consequences of (16. 6p . □ 

We are grateful to Andrea Sportiello for drawing our attention to Lemma 16.11 and 
its proof, and for showing us this elegant proof of Proposition 16.31 

An analogue of Proposition 16.31 holds also for hermitian quaternionic matrices, 
albeit with a different proof: here the determinant is the Moore determinant, and 
"rank" means left row rank (= right column rank); moreover, the polynomial P(x\, . . . , x n ) 
is defined initially by letting xi,...,x n be real numbers. See Proposition IA.81 

6.2 The half-plane property 

Now we take an analytic point of view, so that the commutative ring R will be 
either 1 or C. 

In this subsection we make a slight digression from our main theme, by showing 
that if the are complex hermitian positive-semidefinite matrices (this of course 
includes the special case of real symmetric positive-semidefinite matrices), then the 
determinantal polynomial (16. 7p has the half-plane property. This turns out to be an 
easy extension of the proof of [101 Theorem 8.1(a)]. Indeed, we can go farther, by first 
stating the result in a clean abstract way, and then deducing the half-plane property 
for (16.71) as an immediate corollary. 

Proposition 6.4 Let V be the real vector space Herm(m, C) of complex hermitian 
mxm matrices, and let C C V be the cone Il m (C) of positive-definite matrices. Then 
the polynomial P(A) = detA is nonvanishing on the tube C + iV. 

[Of course, by restriction the same result holds when V is the real vector space 
Sym(m,IR) of real symmetric mxm matrices and C C V is the cone II m (R) of 
positive- definite matrices.] 

Corollary 6.5 Let At, . . . ,A n be complex hermitian positive-semidefinite mxm ma- 
trices. Then the polynomial 



has the half-plane property, that is, either P = or else P(x±, . . . ,x n ) ^ whenever 
ReXj > for all i. 

First Proof of Proposition 16.41 Let A e C + iV, and let ip be a nonzero 
vector in C m . Then the Hermitian form ip*Aip = Y2Tj=i' l t ) iAij'ipj has strictly positive 
real part, and in particular is nonzero; it follows that Aip ^ 0. Since this is true for 
every nonzero ip G C m , we conclude that ker A = {0}, i.e. A is nonsingular; and this 
implies that (and is in fact equivalent to) detA ^ 0. □ 




(6.8) 
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Second Proof of Proposition 16.41 [27J Lemma 4.1]. Write A = P + %Q with P 
positive-definite and Q hermitian. Then P has a positive-definite square root P 1 / 2 , 
and we have 



Proof of Corollary 16.51 If at least one of the matrices Ai,...,A n is strictly 
positive definite, then Y^i=i nes in C + whenever Rex« > for all i. Propo- 
sition [63] then implies that P(xi, . . . , x n ) ^ 0. 

The general case can be reduced to this one by replacing A4 —> A* + el (e > 0) 
and taking e \. 0. By Hurwitz's theorem, the limiting function is either nonvanishing 
on the product of open right half-planes or else is identically zero. □ 

Remarks. 1. The first proof of Proposition 16.41 is an abstraction of the proof of 
Choe et al. [HH Theorem 8.1(a)]; it can also be found in Faraut and Koranyi [5B| 
Lemma X.1.2] and, for V = Sym(m, R), in Hormander [77J p. 85]. Faraut and Koranyi 
furthermore observe that A h-> A^ 1 is an involutive holomorphic automorphism of the 
tube C + iV, having / as its unique fixed point; and these facts hold not only for 
Herm(m, C) but in fact for any simple Euclidean Jordan algebra [56j Theorem X.l.l]. 

2. It is natural to ask whether there is an analogous result for permanents, which 
would extend [101 Theorem 10.2] in the same way that Proposition 16.41 extends [4"0~t 
Theorem 8.1(a)]. The most obvious such extension would be: Let V be the vector 
space W nxn of real m x n matrices (m < n) and let C C V be the cone (0, oo) mxn 
of matrices with strictly positive entries. Then the polynomial Q(A) = per A is 
nonvanishing on the tube C+iV. But this is simply false: for instance, for m = n = 2 



Let us now turn to the main results of this section, namely, Theorems 11.31 and 11.41 
and their consequence, Theorem 11.21 The proof of all these results turns out to be 
surprisingly easy; and all but the "only if" part of Theorems II .31 and II .41 is completely 
elementary. 

In order to make the elementary nature of the proof as transparent as possible, 
we proceed as follows: First we prove the direct ("if") half of Theorem ll.3( a.b) 
by completely elementary methods, without reference to Euclidean Jordan algebras. 
Then we prove Theorem 11.41 we will see that the proof of the direct half of this 
theorem is a straightforward abstraction of the preceding elementary proof in the 
concrete cases; only the converse ("only if") half is really deep. 




to 1. □ 




□ 



6.3 Proof of Theorems D, D and D 
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Proof of the direct half of Theorem I1.3I (a.b). Let us begin with the real 
case and /3 = 1/2. We use the Gaussian integral representation 



(detA)- 1 / 2 = [ exp(-x T Ax) ff (6.10) 
J V n 



where A is a real symmetric positive-definite m x m matrix and we have written 
x = (xi,...,x m ) G M" 1 . This exhibits (detA) -1 / 2 as the Laplace transform of a 
positive measure on II m (R)* = LI m (R), namely, the push-forward of Lebesgue measure 
dx/n m / 2 on M m by the map x i— > xx T . [We remark that this measure is supported 
on positive-semi definite matrices of rank 1.] Alternatively, one can see directly, by 
differentiating under the integral sign in (I6.10p . that the /c-fold directional derivative of 
(det A)' 1 / 2 in directions B\, . . . , Bk G II m (R) has sign ( — l) fc , because each derivative 
brings down a factor — x T 5jX < 0. 

Since a product of completely monotone functions is completely monotone, it fol- 
lows immediately that (det A)~ N ^ 2 is completely monotone for all positive integers 
N. We remark that this can alternatively be seen from the Gaussian integral repre- 
sentation 



r / N N m d (a) 

(detA)~ N / 2 = / exp(-^xM T Ax( Q )) nil^ ( 6 - U ) 

•' a=l a=l j=l V 71 



where we have introduced vectors x^^ G W 71 for a = 1, . . . , N. If we assemble these 
vectors into anmxiV matrix X, then we have exhibited (det A)~ N I 2 as the Laplace 
transform of a positive measure on n m (R), namely, the push- forward of Lebesgue 
measure dX/n mN l 2 on M mx7V by the map X i— > XX T . This measure is supported on 
positive-semidefinite matrices of rank min(N,m). 

Finally, for real values of (3 > (m — l)/2, we use the integral representational 



(det A) 



-P - ^m(m-i)/4 JJr(/?-^) x I e-*^ (det Bf^dB 



_m(m- 
7T v 



B>0 



(6.12) 

where the integration runs over real symmetric positive-definite m x m matrices B, 
with measure dB = f] dBij. This exhibits (detA) -13 as the Laplace transform 

l<i<j<m 

of a positive measure on II m (R). 

The proof in the complex case is completely analogous. For = 1 we use the 
Gaussian integral representation 

(detAr = [eM-- T Az)f[ {dReZ ^ dlmZ ^ (6.13) 



=i 



23 See e.g. [82], [ESI pp. 585-586], [1361 P- 41], [53l Lemma 1], or Theorem 7.2.2 and Corol- 
lary 7.2.4]. 
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where A is a complex hermitian positive-definite matrix, z = (z\, . . . , z m ) G C m and 
denotes complex conjugation. For real values of /3 > m — 1, we use the integral 
representation 24 ! 



m-l "|~ J 

(detA)" " 



vr— 1)/2 ri F (/3-J: 



r m(m- 

j=0 



x y e" tr(AB) (detS^dB , (6.14) 



B>0 



where the integration runs over complex hermitian positive-definite m x m matrices 
B, with measure 

m 

dB = JJdBa 1 [ (dReB ij )(dlmB ij ) . (6.15) 

i=l l<i<j<m 

□ 



Remarks. 1. By applying (16.1 2p for the case m = 2 to A 



Xl + x 3 x 3 

X 3 X 2 + X 3/ 

we obtain after a bit of algebra an alternate derivation of the formula ( I3.2ip for 
(xiX2 + X1X3 + X2X3)~ 13 , valid for (3 > 1/2. 

2. Similarly, by combining (16. 14[) / (I6.15P for the case m = 2 with (11.51) / (ll.6p . we 
obtain after some algebra an explicit formula for £ , 2 4(x) - ^, valid for > 1: 

3 

(xiX 2 + X1X3 + £1X4 + 2:2X3 + x 2 x 4 + x 3 x 4 ) _/3 — v 



00 00 00 00 



2vrr(/3-l)r(/3) 

J dt x J dt 2 J dt 3 J dU e -t^-^2-t 3 x 3 -t iX4 Q(t^t 2 MMf+ 2 (6.16) 
0000 

where 

Q{t u t 2 ,hM) = (hh + hts + hU + t^ + hU + tsU) - (tl + tl + tj + tl) . (6.17) 

This formula provides an explicit elementary proof of the direct half of Corollary 11.61 
- and in particular solves the Lewy-Askey problem — in the same way that (I3.2ip 
provides an explicit elementary proof of the direct half of Proposition 13.81 Note also 
that by setting x 4 = in (16 . 16[) and performing the integral over t 4 , we obtain (I3.2ip . 
Finally, see Corollary 17.81 for a generalization from E 23 and E 24 to E 2n . □ 



As preparation for the proof of Theorem 11.41 let us review some facts from the 
theory of harmonic analysis on Euclidean Jordan algebras (see (56] for definitions and 
further background). 

Let V be a simple Euclidean (real) Jordan algebra of dimension n and rank r, with 
Peirce subspaces Vij of dimension d; recall that n = r + |r(r — 1). It is illuminating 
(though not logically necessary) to know that there are precisely five cases [561 P- 97]: 

24 See e.g. [69l eqns. (5.13)-(5.40)]. Please note that [7J Problem 7.9] and [70l eq. (1.2)] appear 
to contain errors in the normalization constant. 
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(a) V = Sym(m, K), the space ofmxm real symmetric matrices (d = 1, r = m); 

(b) K = Herm(m,C), the space of m x m complex hermitian matrices (d = 2, 
r = m); 

(c) V = Herm(m, H), the space ofmxm quaternionic hermitian matrices (d = 4, 
r = m); 

(d) V = Herm(3, O), the space of 3 x 3 octonionic hermitian matrices (d — 8, r = 3); 
and 

(e) V" = R x M™" 1 (d = 7i - 2, r = 2). 

We denote by (x\y) = tr(xy) the inner product on V, where tr is the Jordan trace 
and xy is the Jordan product @ 

Let Q C V be the positive cone (i.e. the interior of the set of squares in V, or 
equivalently the set of invertible squares in V); it is open, convex and self-dual@ We 
denote by A(x) the Jordan determinant on V: it is a homogeneous polynomial of 
degree r on V, which is strictly positive on Q and vanishes on <9fi@ We have the 
following fundamental Laplace-transform formula \5E\ Corollary VII. 1.3]: for y G Q 
and Re a > (r - l)f = 2 - 1, 

J e- {xly) A(x) a -^ dx = T n (a)A(y)- a (6.18) 
h 



25 In cases (a) and (b), we have (A\B) = tr(AB); in case (c) we have (A\B) = \ tr(AB + BA); 
and in case (e) we have ((x , x)|(y , y)) = 2(x y + x- y). 

26 In cases (a)-(d) the positive cone f2 is the cone of positive-definite matrices; in case (e) it is the 
Lorentz cone {(xq,x): xq > Vx 2 }. 

27 In cases (a) and (b), the Jordan determinant is the ordinary determinant; in case (c) it is the 
Moore determinant (see Appendix in case (d) it is the Freudenthal determinant [5T1H51 11711102) : 
and in case (e) it is the Lorentz quadratic form A(xq,x) — Xq — x 2 . 
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whercfl ^ 

r n (a) = (2tt)^)/ 2 Jjr(a-^) . (6.19) 

Thus, for Re a > (r - 1)| , the function A(x) a -*/r n ( a) is locally integrable on Q 
and polynomially bounded; it therefore defines a tempered distribution TZ a on by 
the usual formula 

TZJip) = — — - I tpfx) A(x) a -r dx for up e S(V) . (6.20) 

n 

Using f)6.18p . a beautiful argument — which is a special case of I. Bernstein's general 
method for analytically continuing distributions of the form Vq [2"2"| |2"3"] — shows that 
the distributions lZ a can be analytically continued to the whole complex a-plane: 

Theorem 6.6 [56, Theorem VII. 2. 2 et seq.] The distributions lZ a can be 
analytically continued to the whole complex a-plane as a tempered- distribution-valued 
entire function of a. The distributions lZ a have support contained in Q and have the 
following properties: 

TZ = 5 (6.21a) 
n a *lZ p = TZ a+p (6.21b) 
A(d/dx)TZ a = n a ^ (6.21c) 

A( X )n a = (n{ a -4)) na+i (6 - 2id) 

(here S denotes the Dirac measure at 0). Finally, the Laplace transform of1Z a is 

(Ol a )(y) = A(y)- a (6.22) 



2<s 



Here dx is Lebesgue measure on the Euclidean space V with inner product ( • | • ) . Thus in 



case (a) we have (B\B) = Y B 2 U + 2 Y B ij> so that 



I J 1 

i—l l<i<j<m 



dx = Y[dB u Y[ ^ dB v = 2 m(m - 1)/4 dB, 

i— 1 l<i<j<m 

m 

showing that (|6. / (|6.19[) agrees with (|6.12[) . Likewise, in case (b) we have (B\B) = Y B u 

i=l 

2 Y \ B ij\ 2 , so that 

l<i<j<m 

m 

dx = Y[dBu Y[ (V^dReB^iV^dlmBij) = 2 m ( m - 1)/2 dB , 

i—l l<i<j<.m 

showing that (f57T5|> / ^HUHD agrees with ([6TT1) . 
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for y in the complex tube Q + ivffi 

The distributions {TZ a }a£C constructed in Theorem 16.61 are called the Riesz distribu- 
tions on the Euclidean Jordan algebra V. 

It is fairly easy to find a sufficient condition for a Riesz distribution to be a 
positive measure: 

Proposition 6.7 J56J, Proposition VII.2.3] (see also [73], Section 3.2] [551125] ) 

(a) For a = fe| with k = 0, 1, ... ,r — 1, the Riesz distribution lZ a is a positive 
measure that is supported on the set of elements of Q of rank exactly k (which 
is a subset of dVt). 

(b) For a > (r — 1)|, the Riesz distribution lZ a is a positive measure that is sup- 
ported on fl and given there by a density (with respect to Lebesgue measure) that 
lies in Lj oc (fi). 

Indeed, part (b) is immediate from the definition ( 16.201) . while part (a) follows by 
reasoning that abstracts the constructions given in (16.101) / (16.111) an d (16.131) above for 
the special cases of real symmetric and complex hermitian matrices J 30 l 

It is a highly nontrivial fact that the converse of Proposition 16.71 also holds: 

Theorem 6.8 [56, Theorem VII. 3.1] The Riesz distribution lZ a is a positive 
measure if and only if a — 0, |, . . . , (r — 1)| or a > (r — 1)|. 

This fundamental fact was first proven by Gindkin [57] (see also [HS]) and is generally 
considered to be deep. However, there now exist two elementary proofs: one that is 
a fairly simple but clever application of Theorem 16.61 and Proposition 16.71 [127] 134"] 
|131[ Appendix], and another that analyzes the integrability of A(x) a ~~ near dQ and 
characterizes those a G C for which 7Z a is a locally finite complex measure [131] . 

Using Proposition 16.71 for the case d — 4, we can prove the direct half of Theo- 
rem [L3lc) [i.e. the quaternionic case] by complete analogy with the elementary proofs 
given above for the real and complex cases. Moreover, with Theorem 16.81 in hand, the 
proof of Theorem 11.41 (and hence also of the converse half of Theorem I1.3P becomes 
utterly trivial: 

29 The property (|6.21d[) is not explicitly stated in [56], but for Re a > (r — 1)| it is an immediate 
consequence of (|6.19p / (|6.20[) . and then for other values of a it follows by analytic continuation (see 
also [721 Proposition 3.1 (hi) and Remark 3.2]). 

30 Thus, for integer N > in the real symmetric (resp. complex hermitian) case, the positive 
measure Rm/2 is supported on the positive-semidefinite matrices of rank min(iV, m) [561 Propo- 
sition VII.2.3] and is nothing other than the push-forward of Lebesgue measure on M mxAr (resp. 
C mxJV ) by the map X H> XX T (resp. X h-> XX*), as discussed above during the proof of the direct 
half of Theorem Ota). For N > m this is a straightforward calculation [13! HI E3 EH [Jj, which 
shows the equivalence of (|6. 1 1|) and (|6.12p [or the corresponding formulae in the complex case]; for 
< N < m — lit follows by comparing (|6. 1 1|) with (|6.22p and invoking the uniqueness of inverse 
Laplace transforms on the space of distributions |126[ p. 306]: see [56l Proposition VII. 2. 4]. 



44 



Proof of Theorem 11.41 Equation (I6.22p writes A(y)~ a as the Laplace transform 
of the distribution lZ a , which is supported on Q (this is the closed dual cone of fl 
since Q is self-dual). By the Bernstein-Hausdorff-Widder-Choquet theorem (Theo- 
rem [22]), it follows that the map y h-> A(y)~ a is completely monotone on Q if and 
only if TZ a is a positive measure; moreover, if lZ a is not a positive measure, then this 
map is not completely monotone on any nonempty open convex subcone Q' C Q. So 
Theorem 11.41 is an immediate consequence of Theorem 16.81 □ 

Remark. The formulae in this section arise in multivariate statistics in con- 
nection with the Wishart distribution [1061 [7] ; in recent decades some statisticians 
have introduced the formalism of Euclidean Jordan algebras as a unifying device 
[23 EH ESI ESI ES]- These formulae also arise in quantum field theory in studying 
the analytic continuation of Feynman integrals to "complex space-time dimension" 
[1321 [2H [53]. □ 

Theorem 11.31 is an immediate consequence of Theorem ll.4| once we use the fact 
that the Jordan determinant coincides with the ordinary determinant on Sym(m, R) 
and Herm(m, C) and with the Moore determinant (see Appendix [A} on Herm(m, HI). 
Theorem 11.21 is in turn an easy consequence of Theorem 11.31 

Proof of Theorem 11.21 If Ax, . . . , A n are strictly positive definite, Theorem Ol 
is an immediate corollary of Theorem 11.31 and the definition of complete monotonicity 
on cones (Definition I2.ip . The general case is reduced to this one by replacing Aj — > 
Ai + el (e > 0) and taking e | 0. □ 

Conversely, if (3 does not lie in the set described in Theorem II. 3[ then the map 
A i — y (det A)~P is not completely monotone on any nonempty open convex subcone of 
the cone of positive-definite matrices. This means, in particular, that if the matrices 
Ax, ... , A n span Sym(m, K) or Herm(m, C) [so that the convex cone they generate 
is open], then the determinantal polynomial ( 16. 8 j) does not have P' 13 completely 
monotone on (0, oo) n . The spanning-tree polynomial of the complete graph K m+ x 
is an example of this situation; and there are two other cases arising from complex- 
unimodular matroids. The following result thus provides a (very) partial converse 
to Corollary II. 8| where part (a) concerns regular [= real-unimodular] matroids and 
part (b) concerns complex-unimodular matroids: 

Corollary 6.9 Let M be a matroid on the ground set E, and let Bm(x) be its basis 
generating polynomial. 

(a) If M = M(K p ) [the graphic matroid of the complete graph K p ] and (3 ^ 
{0, |, 1, |, . . .} U [{p — 2)/2, oo), then B^f is not completely monotone on any 
nonempty open convex subcone of (0, oo) E . 
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(b r= 2) If M = C/2,4 [the uniform matroid of rank 2 on four elements] and (3 {0} U 
[l,oo), then Bjf is not completely monotone on any nonempty open convex 
subcone of (0, 00) E . 

(K=-i) If M = AG (2, 3) [the ternary affine plane] and f3 {0, 1} U [2, 00), then B^f is 
not completely monotone on any nonempty open convex subcone of (0, 00) E . 

Proof. When M is a real-unimodular (resp. complex-unimodular) matroid of rank r 
with n elements, we let B be an r x n real-unimodular (resp. complex-unimodular) 
matrix that represents M; we then define A{ (1 < i < n) to be the outer product 
of the ith column of B with its complex conjugate. We shall show that in the cases 
enumerated above, we can choose B so that the matrices A\, . . . , A n span Sym(r, R) 
[resp. Herm(r, C)]. The result then follows from Theorem 11.31 and the observation 
made immediately after it. 

(a) Let Hp be the directed vertex-edge incidence matrix for K p with an arbitrarily 
chosen orientation of the edges; it is of size p x and is real-unimodular. Then 
K p is represented over R by the matrix H' obtained from H p by deleting one of the 
rows. But we can reorder the columns of H' p to get H' p ' = (I p _i\H p -i) where I p -\ 
is the (p — 1) x (p — 1) identity matrix (and H p _i is defined using the orientation 
of Kp-x inherited from K p ). The corresponding matrices A\, . . . ,A^py obtained by 

taking outer products of the columns of H' p ' with their complex conjugates, are easily 
seen to span Sym(p — 1, R). 

(b r= 2) The matrices A\, . . . , A4 defined in (11 .6p are easily seen to span Herm(2, C). 

(b r=3 ) Write u = e ±27ri/3 ; then the matrix [HSJ p. 597] 

/ 1 1 1 1 1 1 \ 
5= 010110 1 + 0; 1 1 + u \ (6.23) 
\0 1 1 u to l + u 1 + u J 

is easily seen to be complex-unimodular and to represent AG (2, 3). A tedious compu- 
tation (or an easy one using Mathematica or Maple) now shows that the matrices 
Ai, . . . , A 9 are linearly independent, hence span the 9-dimensional space Herm(3, C). 
□ 

Let us remark that the cases enumerated in Corollary l6.9l exhaust the list of regular 
or complex-unimodular matroids (of rank r > 2) for which the matrices A\, . . . , A n 
span Sym(r, R) or Herm(r, C), respectively. Indeed, it is known that a simple rank-r 
matroid that is regular (or, more generally, is binary with no F7 minor) can have at 
most r(r + l)/2 elements; furthermore, the unique matroid attaining this bound is 
M(K r+ i) |110[ Proposition 14.10.3]. See also [16] for an intriguing proof that uses 
the matrices A±, . . . , A n (but over GF{2) rather than C). Likewise, it is known |1121 
Theorem 2.1] that a simple rank-r matroid that is complex-unimodular can have at 
most (r 2 + 3r — 2)/2 elements if r 7^ 3, or 9 elements if r = 3; furthermore, the 
unique matroid attaining this bound is T r (defined in |112j ) when r 7^ 3, or AG(2, 3) 
when r = 3. The only cases in which this size reaches dimHerm(r, C) = r 2 are thus 
r = 1, 2, 3, yielding T\ = Ui t i, T 2 = ?7 2 ,4 and AG (2, 3), respectively. 
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6.4 Combining the determinantal method with constructions 

Let us now combine the ab initio results from the determinantal method (The- 
orems ll.2til.4l and their corollaries) with the constructions from Sections |3] and H] 
(deletion, contraction, direct sum, parallel connection and series connection). For 
graphs we have: 

Proof of Proposition 11.121 For p = 2 the result is trivial, as the graphs 
concerned are precisely those that can be obtained from a forest by parallel extensions 
of edges and additions of loops (Theorem 15. 3j) . For p > 3 (hence (3 > 1/2), the result 
is an immediate consequence of Corollary 11.51 and Propositions 14.2} 14.31 and 14.41 □ 

The matroid generalization of Proposition 11.121 is: 

Proposition 6.10 Fix r > 1, and let M be any matroid (on the ground set E) that 
can be obtained from regular matroids of rank at most r by parallel connection, series 
connection, direct sum, deletion and contraction. Then B A f is completely monotone 
on (0, oo) E for /3 = 0, |, 1, §, . . . and for all real > (r - l)/2. 

The proof is essentially identical to the previous one, but uses Corollary 11.8( a) in 
place of Corollary 11.51 and Propositions I4.6ti4.8l in place of Propositions I4.2ti4.4l 
And we also have: 

Proposition 6.11 Fix r > 1, and let M be any matroid (on the ground set E) 
that can be obtained from regular matroids of rank at most 2r — 1 and complex- 
unimodular matroids of rank at most r by parallel connection, series connection, direct 
sum, deletion and contraction. Then B]f is completely monotone on (0, oo) E for 
j3 = 0, 1, 2, 3, . . . and for all real > r — 1. 

The proof is again identical, but uses both parts of Corollary 11.81 instead of only 
part (a). 

Propositions 11.121 16.101 and 16.111 give a rather abstract characterization of the 
class of graphs or matroids that they handle, and so it is of interest to give a more 
explicit characterization. Let us start with the graph case. We say that a graph G is 
minimally 3-connected if G is 3-connected but, for all edges e of G, the graph G\e 
is not 3-connected. We then have: 

Proposition 6.12 Let Q p be the class of graphs obtained from K p by disjoint union, 
gluing at a cut vertex, series and parallel connection, deletion and contraction of 
edges, and deletion of vertices. Then, for p > 3, Q p is minor-closed, and the minimal 
excluded minors for Q p are the minimally 3-connected graphs on p + 1 vertices. 

Proof. Let l-i p be the class of graphs with no minimally 3-connected minor on p+ 1 
vertices. It is clear that Q p is minor-closed, so our aim is prove that Q p = ~K P . 
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We first show that Q p C H p . It is clear that 7i p is minor-closed, and is closed under 
disjoint union ("0-sum") and gluing at a cut vertex ("1-sum"). It is easily checked 
that if a graph G is obtained by parallel connection of graphs Gi and G2, then any 
3-connected minor of G is a minor of either G\ or G 2 ; it follows that H p is closed 
under parallel connection. Since series connection can be obtained by combining the 
other operations (exploiting K 3 G Hp), we conclude that % p is also closed under series 
connection. Finally, we note that K p G fi p , and as Q p is the closure of {K p } under 
these operations, we see that Q p C % p . 

We now show that % p C Q p . For suppose otherwise, and choose G G 7i p \ Q p 
with the minimal number of vertices. Clearly G is 2-connected and has at least 
p + 1 vertices. If G is not 3-connected, then G has a cutset {x, y} and there is a 
decomposition G = G\ U G2 where G\ and G2 are connected graphs with at least 3 
vertices and V(G\) fl V(G 2 ) = {x, y}. Now let G[, G' 2 be the graphs obtained from 
G\, G 2 by adding the edge xy if not present. Then G[ and G' 2 are both minors of G 
(obtained by contracting the other side to a single edge). As % p is minor-closed, we 
have G^G'2 G Tip. Therefore, by minimality of G we have GJ, G' 2 G Q p . But G can 
be obtained by taking a parallel connection of G' x and G' 2 along and then deleting 
the edge xy if necessary, yielding G G Q p , contrary to hypothesis. We conclude that 
G must be 3-connected. 

We now use the fact that every 3-connected graph other than K4 has an edge 
that can be contracted to produce another 3-connected graph |l8j Lemma 3.2.4]. 
Contracting suitable edges, we see that G has a 3-connected minor on p + 1 vertices, 
which in turn (deleting edges if necessary) contains a minimally 3-connected minor 
on p + 1 vertices. But this contradicts G G H p . □ 

The matroid case is analogous. We refer to |110[ Chapter 8] for the definitions of 
3-connectedness and minimal 3-connectedness for matroids. The following result and 
its proof are due to Oxley 

Proposition 6.13 Let T be a class of matroids that is closed under minors, direct 
sums and 2-sums. Let r > 1, let T T = {M G T: rank(M) < r}, and let T' T denote 
the class of matroids that can be obtained from T r by direct sums and 2-sums. Then 
T' T is closed under minors (and of course also under direct sums and 2-sums), and it 
consists of all matroids in J 7 that have no minor, where 

= {M G T: M is minimally 3-connected and rank(M) = r + 1} . (6.24) 

Proof. Consider first the case r — 1. The unique minimally 3-connected matroid of 
rank 2 is £/2,3- If ^2,3 G J 7 , then the result clearly holds; if f/2,3 ^ J 7 , then J r [ = T 
and the result again holds. 

Now assume r > 2, and let "H r denote the class of matroids in T that have no 
J-'r+i minor. Clearly "H r is minor-closed, and our goal is to show that T' r = H r - 

We first show that T' r C % r . It is easy to see |110[ Propositions 4.2.20 and 8.3.5] 
that if a matroid M is either a direct sum or a 2-sum of matroids M\ and M 2 , then 
any 3-connected minor of M is isomorphic to a minor of either M\ or M 2 ; it follows 
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that H r is closed under direct sum and 2-sum. Since T r C % r , and T' r is the closure 
of T T under direct sum and 2-sum, we see that T' r ^Tir- 

We now show that ^ r C J^. For suppose otherwise, and choose MeH r \JJ with 
the minimal number of elements. It is not hard to see that M must be 3-connectedJ^I 
Note also that rank(M) > r since M E 7 \ T' r . 

We now use Tutte's Wheels and Whirls Theorem [1101 Theorem 8.8.4], which says 
that every 3-connected matroid N that is not a wheel or a whirl has an element e 
such that either iV \ e or N/e is 3-connected (or both). On the other hand, if N 
is a wheel or a whirl, then by contracting a rim element and deleting one of the 
spokes adjacent to that rim element, we obtain a 3-connected minor N' of N such 
that rank(jV') = rank(iV) - 1. 

So we apply this argument repeatedly to M until we obtain a 3-connected minor 
M' of M with rank(M') = r + 1. We then delete elements from M' while maintaining 
3-connectedness until we arrive at a minimally 3-connected matroid M". Therefore 
M has a minor M" G contradicting the hypothesis that M 6 T-L r . □ 

Proposition 16.131 with J 7 = regular matroids gives an excluded-minor characteri- 
zation of the class of matroids handled by Proposition 16. 101 We leave it as a problem 
for readers more expert in matroid theory than ourselves to provide an analogous 
characterization for Proposition 16.111 



7 Quadratic forms: Proof of Theorem 11.9 



In this section we consider quadratic forms (= homogeneous polynomials of de- 
gree 2). We begin by proving an abstract theorem giving a necessary and sufficient 
condition for such a quadratic form to be nonvanishing in a complex tube C + iV; in 
the special case C = (0, oo) n this corresponds to the half-plane property. Finally, we 
prove Theorem 11.91 



7.1 The half-plane property 

In this subsection we proceed in three steps. First we study the analytic geometry 
associated to a symmetric bilinear form B on a finite-dimensional real vector space V 
(Lemma 17. II) . Next we extend the quadratic form Q(x) = B(x, x) to the complexified 
space V + iV and study the values it takes (Proposition 17. 2p . Finally we introduce the 
additional structure of an open convex cone C C V on which Q is assumed strictly 
positive (Corollary 17.31 and Theorem I7.4p . 

So let V be a finite-dimensional real vector space, and let B: V x V — > R be a 
symmetric bilinear form having inertia (n+,n_,no). Define <S+ = {x: B(x,x) > 0} 



31 If M were a direct sum or a 2-sum of matroids M\ and M2, each having at least one element, 
then Mi and M2 would be minors of M [1101 Proposition 7.1.21], hence M\^Mi G W r because H r 
is minor-closed. But Mi and M2 cannot both belong to T' r because T' r is closed under direct sum 
and 2-sum and M ^ T' r . Therefore Mi or M2 would be a counterexample to the minimality of M. 
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and «S_ = {x: B(x, x) < 0}. Clearly S + and <S_ are open cones (not in general convex 
or even connected). Indeed, S + and 5_ are never convex (except when they are 
empty) because x G S± implies — x G S± but manifestly ^ S±. 

Many of our proofs will involve choosing a basis in V (and hence identifying V 
with MJ 1 ) in such a way that B takes the form 

B (x,y) = ^xiVi - ^2 x iVi ■ ( 7A ) 

i=l i=n_|_+l 

Moreover, whenever 5 + / we can choose the basis in this construction such that 
the first coordinate direction lies along any desired vector x £ S + : that this can be 
done follows from the standard Gram-Schmidt proof of the canonical form (17. ip . 
Elementary analytic geometry gives: 

Lemma 7.1 Let V be a finite- dimensional real vector space, and let B be a symmetric 
bilinear form on V having inertia (n + , n_, n^). 

(a) Ifn + = 0, then S+ = 0. 

(b) Ifn + = 1, then S + is a nonempty disconnected open non-convex cone, and there 
exists a nonempty open convex cone C such that S + = CU (— C) and Cfl (— C) = 
0. [Clearly C is uniquely determined modulo a sign.] Moreover, B(x,y) > 
B(x, x) 1 ^ 2 B(y,y) 1 ^' 2 > whenever x,y G C; and B(x,y) 2 > B(x,x)B(y,y) 
whenever x G S + and y G V (or vice versa). 

(c) Ifn + >2, then S + is a nonempty connected open non-convex cone. Moreover, 
for each x G S + , the set 

%{x) = {y G V: B(x, y) 2 < B(x, x) B(y, y)} (7.2) 

is a nonempty open non-convex cone that is contained in S + and has a nonempty 
intersection with every neighborhood of x; moreover we can write 

7+(x) = {y EV: IRx+Ry is a two-dimensional subspace contained within 5 + U{0} } 

(7-3) 

[that is, it is a two-dimensional subspace on which B is positive- definite]. 
Analogous statements hold for 5_ when = 0, n_ = 1 or n_ > 2. 

Proof, (a) is trivial. 

(b) Assume that B takes the canonical form (17.11) with n + = 1, and define C to 
be the "forward light cone" 

C = {x G W 1 : x\ — x\ — . . . — x 2 n _ +l > and x x > 0} . (7.4) 

It is immediate that S + = C U (— C) and C PI (— C) = 0, and the statements about S + 
follow easily. 
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Now consider any x,y G C and define 

g(a) = B(x + ay,x + ay) = B(x,x) + 2aB(x,y) + a 2 B(y,y) . (7.5) 

We have g(0) = B(x, x) > 0; but for the special value a* = —x\jy\ the vector x + a±y 
has its first component equal to zero and hence = B(x+a+y, x+a+y) < 0. So the 
quadratic equation g(a) = has a real solution, which implies that its discriminant 
is nonnegative, i.e. that B(x,y) 2 > B(x,x) B(y,y). 

Next assume that x G S + and y G V. If B(x,x) = or B(y,y) < 0, the assertion 
is trivial; so we can assume that x,y G S + . By the replacements x — >■ ±x and y — > ±y 
(which do not affect the desired conclusion) we may assume that x, y G C. But in 
this case the desired inequality has already been proven. 

Finally, using B(x, y) > for x, y G C it is easily checked that C is convex. 

(c) Clearly S + is nonempty; and as explained earlier it is non-convex. To prove 
that S + is connected, we can assume that B takes the form (17. ip with n + > 2. 
It is now sufficient to find a path in S + from an arbitrary vector x G S + to the 
vector e\ = (1,0, ...,0). But this is easy: first move coordinates Xi with % > n + 
monotonically to zero [this increases B(x,x) monotonically and hence stays in «S+]; 
then rotate and scale inside the subspace spanned by the first n + coordinates to 
obtain e\. 

Now assume that B takes the canonical form (17.11) with n + > 2 and with the given 
vector x G S + lying along the first coordinate direction. Then an easy computation 
shows that y = (y lt y 2 , . . . , y n ) belongs to T+(x) if and only if y' = (0, y 2 , . . . , y n ) 
belongs to S + . Therefore, the preceding results (b,c) applied with n + replaced by 
n + — 1 show that 7+(x) is a nonempty open non-convex cone, which is obviously 
contained in S + ; and by taking y' small we see that 7+ (x) meets every neighborhood 
of x. Moreover, + ~Ry = Rx + Ry' is a two-dimensional subspace if and only if 
y' 7^ (i.e. y ^ Rx); and since B(x, y') = 0, we have Ra; + M.y' C S + U {0} if and only 
if y' G 5 + U {0}. □ 

Remark. Note the sharp contrast between parts (b) and (c): in the latter case, 
given any x E S + there is a nonempty open cone of vectors y satisfying the Schwarz 
inequality (strictly) with x\ while in the former case all vectors y G V satisfy the 
reverse Schwarz inequality with x. □ 

We now consider the quadratic form Q(x) = B(x, x), extended to the complexified 
space V + iV in the obvious way: Q(x + iy) = B(x,x) — B(y,y) + 2iB(x,y). We 
want to study the values taken by Q in the complex tubes S + + iV and 5_ + iV, and 
in particular the presence or absence of zeros. We write EI to denote the open right 
half-plane {C G C: ReC > 0}. 

Proposition 7.2 Let V be a finite- dimensional real vector space, let B be a sym- 
metric bilinear form on V having inertia (n + ,n_,n ) , and let Q be the associated 
quadratic form extended to V + iV . 
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(a) If n + = 1, then for every x G S + we have Q[Mx] = Q[S + + iV] = C \ (—00, 0]. 
In particular, Q is nonvanishing on S + + iV . 

(b) If n + > 2, then for every x G S + and y G T+(x) [recall that 7~+{x) is nonempty] 
we have Q[[l,oo)x + i(M.x + Ry)] = C. In particular, for each x G <S+ and 
y G 7+(x) there exists z G + My such that Q(x + iz) =0. 

Analogous statements hold for S- when n_ = 1 or n_ > 2. 

Proof, (a) If x G S + , then Q(C^) = C 2 Qi x ) can take any value in C \ (— oo,0] 
as C ranges over EL On the other hand, if x G <S+ and y G V, then Q(x + zy) = 
x) — B(y, y) + 2iB(x, y) cannot take a value in (—00, 0], as B(x, y) — implies 
by Lemma mT b) that y ^ S + , i.e. B(y,y) < 0. 

(b) Given a; G 5 + and ?/ G 7+ (a;), the vector = y+fix with = —B(x, y)/B(x, x) 
satisfies B(x,y') = and Q(y') > 0. Therefore 

Q{\x + i{ax + I3y')) = (A + ia) 2 Q(x) - /3 2 Q(y') , (7.6) 

and this is easily seen to take all complex values as A ranges over [1, 00), a over M, 
and (3 over (0, 00). □ 

Now we introduce the additional structure of an open convex cone C C V on 
which Q is assumed strictly positive (i.e. C C <S+). The hypotheses in the following 
result are identical to those of Theorem 11.91 

Corollary 7.3 Let V be a finite- dimensional real vector space, let B be a symmetric 
bilinear form on V having inertia (n + ,ri-,no), and define the quadratic form Q(x) = 
B(x,x). Let C C V be a nonempty open convex cone with the property that Q(x) > 
for all x G C. Then n + > 1, and moreover: 

(a) Ifn + = 1, then either C C C orC C — C [where C is defined as in Lemma YTT^ b )], 
and we have B(x, y) > for all x,y G C . 

(b) Ifn + > 2, then for every x G C we have Q[[l, oo)x + iV] = C, and in particular 
there exists z G V such that Q(x + iz) = 0. 

Proof, (a) Let us use the canonical form (17. ip with n + = 1 and C defined by (I7.4p . 
By hypothesis the cone C is contained within S + = C U (— C); but since C is convex 
we must in fact have either C C C or C C — C: for otherwise C would contain a 
point with x\ — and hence B(x,x) < 0@ The remaining statements follow from 
Lemma [7. 1( b) and Proposition 17.2( a). 



32 Alternatively, we can argue (assuming for simplicity that no = 0) that a convex set S contained 
in 5+ must either lie on a single line through the origin (which is obviously impossible if S is open) 
or else be contained in either C or —C. For if we had 1 e SflC and j/eSfl — C with x and y not on 
the same line through the origin, then there would exist A € (0, 1) such that z — Xx + (1 — X)y has 
Z\ = and z ^ 0, hence 2 ^ iS+. 
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(b) Since C C S + , the statements follow from Lemma l7Tl7 c) and Proposition ^. 2( b). 

□ 

We can also summarize our results in a way that extends [4"0l Theorem 5.3] from 
C = (0, oo) n to general open convex cones C: 

Theorem 7.4 Let V be a finite- dimensional real vector space, let B be a symmetric 
bilinear form on V having inertia (n + , n_, n ), and define the quadratic form Q{x) = 
B(x,x). Let C C V be a nonempty open convex cone with the property that Q(x) > 
for all x G C. Then n + > 1, and the following are equivalent: 

(a) n + = 1. 

(b) Q is nonvanishing on the tube C + iV . 

(c) Q[C + tV] = C\(-oo,0]. 

(d) If x,y eV with Q(x) > 0, then B(x, yf > B(x, x) B(y, y) . 

(e) Ifx,ye C, then B[x,yf > B(x,x) B(y,y). 

Proof. This follows immediately by putting together the statements from Lemma I7TT1 
Proposition 17.21 and Corollary 17.31 if n + = 1, then (b,c,d,e) are all true; and if n + > 2, 
then (b,c,d,e) are all false. □ 

Remark. See pEDJ Remark 1 after the proof of Theorem 5.3] for some of the 
history of this result in the traditional case C = (0, oo) n . □ 

7.2 Proof of Theorem [T791 

Proof of Theorem 11.91 We are concerned with the complete monotonicity of 
Q~P , where Q(x) = B(x, x). As always, it suffices to consider (3 > 0, because complete 
monotonicity trivially holds when = and never holds when /3 < (because Q grows 
at infinity). 

We assume that B takes the canonical form ( 17.1 j) on V = M. n , and we consider 
separately the three cases: 

(a) The case n+ = 1, n_ = is trivial: we have Q(x) = x\, and the convex 
cone C must be contained in one of the half-spaces {x\ > 0} or {x\ < 0}. The map 
x i— > Q(x)^ 13 = x^~ 2/3 is clearly completely monotone on each of these two half-spaces. 

(b) Next consider the case n + = 1, rt_ > 1: here ( 17. ip is the Lorentz form in one 
"timelike" variable and n_ "spacelike" variables, and we have S + = CU (-C) where C 
is the forward light cone ( I7.4p . By Corollary 17.3( a) we have either C C C or C C —C; 
let us suppose the former. 

Let us now show that if > (n_ — l)/2, then the map x h-> Q(x)~^ is completely 
monotone on C. The variables x n _ + 2, . . . ,x n play no role in this, so we can assume 
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without loss of generality that no = 0, i.e. n = n_ + l. For (3 > (n_ — 1)/2 = (n— 2)/2, 
the desired complete monotonicity then follows from the integral representation |121[ 
pp. 31-34] [SB eqns. (24)-(28)] 



/ z z z \ 

(a^ — x 2 — ... — x n j 



2\-/3 



7r (r l -2)/2 2 2^1 r( ^ r ^_n_^ 



1-1 

X 



C 

valid for /3 > (n — 2)/2, which explicitly represents Q{x)~ 13 as the Laplace transform 
of a positive measure supported on the closed forward light cone C* = C. For (3 = 
(n — 2)/2 the result follows by taking limits. 

Conversely, let us show that if (3 ^ {0} U [(n_ — l)/2,oo), then the map x i— )■ 
is not completely monotone on any nonempty open convex subcone C" C C. 
For suppose that this map is completely monotone on C for some such j3: then 
by the Bernstein-Hausdorff-Widder-Choquet theorem ( Theorem 12 . 2 p . we must have 
(assuming again without loss of generality that n = Ojfj 

(x\-x\- x 2 n )- p = J e -(*iw-*2W-.— *»w0 dfjip(y) (7.8) 

for some positive measure /i^ supported on (C)*. Now, any such measure must 
clearly be Lorentz- invariant and homogeneous of degree 2/3 — n (this follows from the 
uniqueness of inverse Laplace transforms). Furthermore, must be supported on C, 
for otherwise the support would contain a spacelike hyperboloid {y G W 1 : yf — y^ — 
. . . — y\ = A} for some A < 0, whose convex hull is all of ~&L n and hence not contained 
in the proper cone (C)*. On the other hand, every Lorentz-invariant locally-finite 
positive measure on M n that is supported on C is of the form |120[ Theorem IX. 33, 
pp. 70-76] 

jj, = c5q + / dVL m dp(m) (7.9) 



m>0 



where c > 0, 5 denotes the point mass at the origin, p is a positive measure on [0, oo), 
and dQ m = 5(y 2 — m 2 ) dy is the unique (up to a constant multiple) Lorentz-invariant 
measure on the "mass hyperboloid" 

H m = {yeW'-.yl-yl- . .. -y 2 n = m 2 and yi > 0} . (7.10) 

(When n = 2, we consider m > only, as there is no locally-finite Lorentz-invariant 
measure on M 2 that is supported on Hq.) A measure u of the form ( 17. 9p is homoge- 
neous in precisely three cases: 



33 If the map x (-> Q(x)~@ is completely monotone on C , then it is also completely monotone on 
the cone C" obtained by projecting C onto the first n + +n_ coordinates (since Q(x) is independent 
of the last n coordinates). 
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(a) c = 0, dp(m) = const x m x 1 dm with A > 0: here \x is homogeneous of degree 
A - 2. 

(b) c = 0, p = const x 5 [for n > 2 only]: here /i is homogeneous of degree —2. 

(c) c > 0, p = 0: here is homogeneous of degree — n. 

This proves that a positive measure /ig can exist on/?/ if (3 = or (3 > (n — 2)/2. 

(c) Finally, consider the case n + > 1. By Corollary 17.3( b). Q has zeros in the 
tube C + iV for every nonempty open convex subcone C C 5 + (indeed, for every 
nonempty subset C C «S+). We conclude by Corollary 12.31 that cannot be 

completely monotone on C" for any > 0. □ 

Remark. More can be said about the integral representation (17. 7p . It turns out 
that the quantity 



n 



R B = ^ n -^^- 1 T{p)Y{p-~ 



{yl-yt-----ylYMyzC\, (7.nj 



which is initially defined for f3 > (n — 2)/2 as a positive measure [or for Re /3 > 
(n — 2)/2 as a complex measure] on R™ (and which is of course supported on C), 
can be analytically continued as a tempered-distribution-valued entire function of 
[5Tj [551 Theorem VII. 2. 2]: this is the Riesz distribution Tip on the Euclidean Jordan 
algebra R x R™ _1 @ The integral representation 



^1 -^2 • • • x nJ 



y e -(*i»i-* a «2-...-«„w n ) dy ( 7 _ 12 ) 



where x G C then holds for all complex (3, by analytic continuation. However, the 
distribution Hp is a positive measure if and only if either (3 = or (3 > (n — 2)/2. 
This follows from general results of harmonic analysis on Euclidean Jordan algebras 
(i.e. Theorem 16. 8p . but we have given here a direct elementary proof. □ 



7.3 Explicit Laplace-transform formula for inverse powers of 
a quadratic form 

By a simple change of variables, we can generalize (17. 7p to replace the Lorentz 
matrix L n = diag(l, — 1, . . . , — 1) by an arbitrary real symmetric matrix A of iner- 
tia (n + ,n_,n ) = (l,n — 1,0). This will provide, among other things, an explicit 
Laplace-transform formula for S 2 ,n( x ) - ^ that generalizes the formulae f)3.2ip and 
(I6.16P obtained previously for n = 3,4, respectively. 

34 A slightly different normalization is used in [56j . arising from the fact that the Jordan inner 
product on R x M™ -1 is ((xo,x)|(j/o,y)) = 2(xoVo + x • v) : this has the consequence that dxjordan = 
2™/ 2 (inordinary, and also the Laplace transform is written with an extra factor 2 in the exponential. 
The change of sign from xoyo — x • y to xoyo + x • y is irrelevant, because the Riesz distribution 
IZfj(y) is invariant under the reflections yi H > —yi for 2 < i < n. 
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So let A be a real symmetric n x n matrix with one positive eigenvalue, n — 1 
negative eigenvalues, and no zero eigenvalues. We first need a slight refinement of 
Lemma [7.1( b) to take advantage of the fact that we now have n = 0; for simplicity 
we state it in the "concrete" situation V = M. n . 

Lemma 7.5 Fix n > 2, and let A be a real symmetric nxn matrix with one positive 
eigenvalue, n — 1 negative eigenvalues, and no zero eigenvalues. Then there exists a 
nonempty open convex cone C C M. n (which is uniquely determined modulo a sign) 
such that 



Cn(-C) 
Cn(-C) 
c n {-C) 
C n (-O 

{y: y T Ay > 0} 
{x: x T A" 1 x > 0} 



where 

is the open dual cone 



C = {x: x T y > for 
to C. 




- {0} 

= 

= {0} 

= CU(-C) 

ally GC\{0}} 



(7.13a) 
(7.13b) 
(7.13c) 
(7.13d) 
(7.13e) 
(7.13f) 

(7.14) 



Proof. We can write A = S T L n S where S is a nonsingular real matrix. Then the 
claims follow easily from the corresponding properties of the Lorentz quadratic form. 
□ 



Proposition 7.6 Fix n > 2, and let A be a real symmetric nxn matrix with one 
positive eigenvalue, n — 1 negative eigenvalues, and no zero eigenvalues; and let C be 
the open convex cone from Lemma \7. 5\ Then for > (n — 2)/2 we have 



(^A- 1 *) 



n-2 



{y'Ay f-2 dy 



det^l 1 / 2 



x 



(7.15) 



for x G C 



Proof. Write A = S T L n S where S is a nonsingular real matrix, and make the 
changes of variable y = Sy' and x = L n S~ T x.' in ( 17. 7ft . Then dy = (det^l^y' 
where | det S\ = | det^4| 1//2 ; and the formula ( 17. 15ft follows immediately after dropping 
primes. □ 
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Let us now specialize to matrices A of the form A = XE n — fil n , where I n is the 
n x n identity matrix and E n is the n x n matrix with all entries 1. Then A has 
eigenvalues nX — /i, — /i, . . . , — /i, hence has inertia (n+, n_, n Q ) = (1, n — 1, 0) provided 
that /i > and A > /i/n; and in that case we have A' 1 = \'E n — n'I n where 

A' = - 7 — , // = - . (7.16) 

[The map (A,/i) (-> (A',//) is of course involutive.] Furthermore, we have 

y T Ay = 2A£ 2 , n (y) + (A-/x)||y|| 2 (7.17) 

n 

where E 2tn (y) = ^2 UiUj and ||y|| 2 = vli an d analogously for x T yl _1 x. 

l<i<j<n i=l 

Corollary 7.7 Fza; n > 2, > and A > /i/n. Then for j3 > (n — 2)/2 we have 
2A „ , , (n-l)A- M|| ^ ||2 ^ ; 

-fr2,n(Xj 



/z(nA — /x) ' /j(n\ — /i) 



n-2- n 1 



^"-^(nA-zx) 1 / 2 x 



| e -* T * (2AE 2 , n (y) + (A - ^\\y\\ 2 )^ 2 dy (7.18) 



/or xgC 



Specializing further to the case X — 2/(n — 1), fi — 2 (corresponding to X' — fi' — 
1/2), we obtain: 

Corollary 7.8 (Laplace-transform formula for E^n ) Fix n > 2. Then for (3 > 
(n — 2)/2 we have 

re— 1 _„ n 

= {n ~ l)2 ~ n 2^ / e - xTy K(y) + ^llyllf 

(27r)(- 2 )/ 2 r(/3)r(/3-^) ^ v 2 y 

(7.19) 

forxeC D (0,oo) n . 

For n = 2 this is elementary; for n — 3 it reduces to (13. 21ft ; and for n = 4 it reduces 
to (I6.16p . The formula ( I7.19P provides, in particular, an explicit elementary proof of 
the direct ("if") half of Corollary Offl 



8 Positive-definite functions on cones 

In this section we recall briefly the theory of positive-definite functions (in the 
semigroup sense) on convex cones, which closely parallels the theory of completely 
monotone functions developed in Section [2j We then apply this theory to powers of 
the determinant on a Euclidean Jordan algebra. 



57 



8.1 General theory 

Here we summarize the basic definitions and results from the theory of positive- 
definite functions on convex cones and, more generally, on convex sets. A plethora of 
useful additional information concerning positive-definite (and related) functions on 
semigroups can be found in the monograph by Berg, Christensen and Ressel [2T] . 

Definition 8.1 Let V be a real vector space, and let C be a convex cone in V . Then 
a function f:C — > K is termed positive-definite in the semigroup sense if for all 
n > 1 and all x±, . . . ,x n G C , the matrix {f{x{ + Xj)}™ J=1 is positive-semide finite; or 
in other words, if for all n > 1, all xi, . . . , x n G C and all c±, . . . , c n G C we have 

n 

^clCjfixi + Xj) > 0. (8.1) 

Similarly, a function f:C + iV^-Cis termed positive-definite in the involutive- 
semigroup sense if for all n > 1 and all x%, . . . , x n G C + iV , the matrix 
{f(xi + ~Xj)}ij=i is positive-semide finite. 

Theorem 8.2 Let V be a finite- dimensional real vector space, let C be an open con- 
vex cone in V , and let f: C — > JR.. Then the following are equivalent: 

(a) f is continuous and positive-definite in the semigroup sense. 

(b) f extends to an analytic function on the tube C + iV that is positive- definite in 
the involutive- semigroup sense. 

(c) There exists a positive measure \i on V* satisfying 

f(x) = J e-^d^l) (8.2) 

for all x G C. 

Moreover, in this case the measure /i is unique, and the analytic extension to C + iV 
is given by 

Please note that the completely monotone functions (Theorem 12. 21) correspond 
to the subset of positive-definite functions that are bounded at infinity (in the sense 
that / is bounded on the set x + C for each x G C), or equivalently decreasing (with 
respect to the order induced by the cone C), or equivalently for which the measure 
fi is supported on the closed dual cone C* [rather than on the whole space V* as 
in Theorem 18.2( c)], See |108|. Lemma 1, p. 579] for a direct proof that complete 
monotonicity implies positive-definiteness. 

We remark that the hypothesis of continuity (or at least something weaker, such 
as measurability or local boundedness) in Theorem 18.2( a) is essential, even in the 
simplest case V = K and C = (0, oo). Indeed, using the axiom of choice it can easily 
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be shown [21 pp. 35-36, 39] that there exist discontinuous solutions to the functional 
equation p(x + y) = p(x)p(y) for x, y G (0, oo), and any such function is automatically 
positive-definite in the semigroup sense. However, any such function is necessarily 
non-Lebesgue-measurable and everywhere locally unbounded [21 pp. 34-35, 37-39]. 

Theorem 18.21 is actually a special case of a more general theorem for open convex 
sets that need not be cones. We begin with the relevant definition [681 P- x ] : 

Definition 8.3 Let V be a real vector space. IfC C.V is a convex set, then a function 
f: C — y R is termed positive-definite in the convex-set sense if for all n > 1 and all 
x\,...,x n G C, the matrix {/( X '^ J )}7j=i ^ s positive-semidefinite. More generally, 
if C C V + iV is a conjugation-invariant convex set, then a function f:C^-Cis 
termed positive-definite in the involutive-convex-set sense if for all n > 1 and all 
G C, the matrix {f )}1 j=i ^ positive-semidefinite. 

Note that if C is in fact a convex cone, then a function f:C — > R is positive- 
definite in the convex-set sense if and only if it is positive-definite in the semigroup 
sense. So this concept is a genuine generalization of the preceding one. 

Theorem 8.4 Let V be a finite- dimensional real vector space, let C C. V be an open 
convex set, and let /: C — > R. Then the following are equivalent: 

(a) f is continuous and positive-definite in the convex-set sense. 

(b) f extends to an analytic function on the tube C + iV that is positive- definite in 
the involutive-convex-set sense. 

(c) There exists a positive measure p on V* satisfying 



for all x G C. 

Moreover, in this case the measure p is unique, and the analytic extension to C + iV 
is given by h8.3\) . 

Theorem 18.41 was first proven by Devinatz [46], using the spectral theory of 
commuting unbounded self-adjoint operators on Hilbert space (he gives details for 
dim\^ = 2 but states that the methods work in any finite dimension); see also 
Akhiezer [H pp. 229-231] for the special case in which C is a Cartesian product of 
open intervals. A detailed alternative proof, based on studying positive-definiteness 
on convex sets of rational numbers as an intermediate step [15], has been given by 
Glockner [6"%1 Proposition 18.7 and Theorem 18.8], who also gives generalizations to 
infinite-dimensional spaces V and to operator-valued positive-definite functions. See 
also Shucker |128[ Theorem 4 and Corollary] and Glockner [681 Theorem 18.8] for the 
very interesting extension to convex sets C that are not necessarily open (but have 
nonempty interior): in this latter case the representation (18. 31) does not imply the 




(8.3) 
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continuity of / on C, but only on line segments (or more generally, closed convex 
hulls of finitely many points) within C . But with this modification the equivalence 
(a') <^=> (c) holds. 

Surprisingly, we have been unable to find in the literature any complete proof 
of Theorem 18.21 except as a corollary of the more general Theorem 18.41 But see 
[68} Theorem 16.6] for a version of Theorem 18.21 for the subclass of positive-definite 
functions that are a-bounded with respect to a "tame" absolute value a. 

It would be interesting to try to find simpler proofs of Theorems 18.21 and 18.41 

8.2 Powers of the determinant on a Euclidean Jordan algebra 

We can now deduce analogues of Theorems 11.31 and 11.41 in which complete mono- 
tonicity is replaced by positive-definiteness in the semigroup sense. For brevity we 
state only the abstract result in terms of Euclidean Jordan algebras: 

Theorem 8.5 Let V be a simple Euclidean Jordan algebra of dimension n and rank r, 
with n = r + ^r(r — 1); let f2 C V be the positive cone; let A: V — > R be the Jordan 
determinant; and let (3 e R. Then the following are equivalent: 

(a) The map x i— > A(x) - ^ is positive-definite in the semigroup sense on Q. 

(b) The map x h- > A(x) _/3 is positive- definite in the semigroup sense on some 
nonempty open convex subcone Q' C Q. 

(c) /3 € {0,f,...,(r-l)f}U((r-l)f,oo). 

Theorem 18.51 is an immediate consequence of facts about Riesz distributions - 
namely, the Laplace-transform formula (16.221) and Theorem 16.81 — together with 
Theorem 18.21 Indeed, the proof of Theorem 18.51 is essentially the identical to that of 
Theorem 11.41 but using Theorem 18.21 in place of the Bernstein-Hausdorff-Widder- 
Choquet theorem. The point, quite simply, is that our proof of the "converse" half 
of Theorem 11.41 used only the failure of positivity of the Riesz distribution, not any 
failure to be supported on the closed dual cone (indeed, it is always supported there); 
so it proves Theorem 18.51 as well. 

Example 8.6 Let V be the real vector space Sym(m, R) of real symmetric to x m 
matrices, let II m (R) C V be the cone of positive-definite matrices, and let C C U m (R) 
be the subcone consisting of matrices that are also elementwise strictly positive. It 
follows from Theorem 18.51 that the map A i— > (detA) -|S is positive-definite in the 
semigroup sense on C •<==>- it is positive-definite in the semigroup sense on IT m (R) 
-<=>- f3 G {0, 2, 1, |, . . .} U [(rn — l)/2, oo). This disproves the conjecture of Gurau, 
Magnen and Rivasseau [ZD Section 7, Conjecture 1] that the map A h> (detA)^ 13 
would be positive-definite in the semigroup sense on C for all (3 > 0. □ 
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A The Moore determinant for quaternionic her- 
mitian matrices 



In this appendix we review the definition and properties of the Moore determinant 
for quaternionic hermitian matrices. For the most part we have followed Aslaksen [9] 
and Alesker [5J E] ■ 

The (real) quaternions HI are the associative algebra over K. consisting of objects 
of the form 

q = x l + x\i + x 2 j + x 3 k (x , xi, x 2 , x 3 G R) , (A.l) 
equipped with the multiplication law where 1 is the identity element and 

i 2 = j 2 = k 2 = — 1, ij = —ji = k, jk = —kj = i, ki = —ik = j . (A. 2) 

The conjugate of a quaternion q = x l + x\% + x 2 j + x 3 k is q = x l — X\% — x 2 j — x 3 k. 
Note that conjugation is an involutive antiautomorphism of H, i.e. q = q and WU = vu. 
We also write \q\ 2 = qq — qq — Xq + x\ + x 2 + x\. The complex numbers a + hi can 
be identified with the subalgebra of EI consisting of quaternions al + hi + Oj ; + Ok. 

The hermitian conjugate of a matrix M G H nxn is M* = (M) T . The map M h- > 
M* is an involutive antiautomorphism of the algebra of nxn quaternionic matrices. 
A matrix M is called hermitian if M = M*. We denote by Herm(n, H) the set of 
nxn hermitian quaternionic matrices. 

The map ip: H ->■ C 2x2 defined by 

is an injective *-homomorphism of the algebra of quaternions into the algebra of 2 x 2 
complex matrices. It can equivalently be written as 

^(a + bj) = ^ fora,6GC. (A.4) 

More generally, the map M nxn — > C 2nx2n defined by 

V(A + Bj) = ioiA,BeC nxn (A.5) 

is an injective *-homomorphism of the algebra of nxn quaternionic matrices into the 
algebra of 2n x 2n complex matrices. Its image is 

$(r x ") = {Z G C 2nx2n : JZ = 2.1} (A.6) 

where 

Note that *(M*) = *(M)* = -J^(M) T J; in particular, M is hermitian if and only 
if *(M) is. 
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It is convenient now to define $(M) = Jty(M), i.e. 

*(A + Bj) = (~ B A A n ) fovA,BeC n * n . (A. 



-A -B / 

The image of $ is the same as that of i.e. 

$(F X ") = {Y e C 2nx2n : JY = YJ} . (A.9) 

Note that $(/) = J and $(M*) = -$(M) T . In particular, M is hermitian if and 
only if $(M) is antisymmetric, and we have 

$(Herm(n,H)) = {Y G C 2nx2n : F = -F T and JY = YJ} . (A.10) 

Recall now that on the space Antisyni(2n, C) of 2n x 2n complex antisymmetric 
matrices there is a polynomial pf , called the pfaffian, satisfying (pf A) 2 = det A and 
pf J = 1. Indeed, the pfaffian is uniquely defined by these conditions. Here are some 
of its properties: 

Lemma A.l (Properties of the pfaffian) The map pf: Antisym(2n, C) — > C has 
the following properties: 



(a) pf is a homogeneous polynomial of degree n, with integer coefficients^ 

(b) pf J = 1. 

(c) (pf A) 2 =detA 

(d) pf(AAX T ) = (detX)(pf A) for any 2n x 2n matrix X. 

(e) (minor summation formula for pfaffians | |84L 185^1 ) More generally, we 
have 

pf(XAA T ) = (detX*j)(pf A u ) (A.ll) 

I C [2n] 
|/| = 2m 

/or any 2m x 2n matrix X (m < n). Here X+j denotes the submatrix of X with 
columns I ( and all its rows ), and An denotes the submatrix of A with rows and 
columns I. 

See [1331 EQ E31 [721 El E51 E2] for further information on pfaffians. 



35 More precisely, if A — (ay) £ Antisym(2n, C), then pf A is a homogeneous polynomial with 
integer coefficients in the variables {ay}i<i<j< n . Indeed, all the coefficients lie in {—1,0, 1}. 

36 See also [33l Theorem A. 15 and remarks after it] for an alternative proof using Grassmann- 
Berezin integration. 
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Let us also observe the following useful fact about pfaffians: If B is an arbitrary 
n x n complex matrix, then 



1)1 ( -B T ' 1>[ 



d<A i°j q 1 d<A D 



-£\ / JW 7 
/ J {-I Oj {-B T 

(A.12) 

where the second equality used Lemma lA.l( b.d) and the last equality used row (or 
column) interchanges. 

We now define a "determinant" for hermitian quaternionic matrices only: 

Definition A. 2 (Moore determinant of a hermitian quaternionic matrix) 

The Moore determinant of a matrix M £ Herm(n, H) is defined by 

Mdet(M) = pf($(M)) = pf(J*(M)) . (A.13) 



Moore's original definition |105] was very different from this one, but the two def- 
initions can be proven to be equivalent [H pp. 63-64] |10U[ pp. 141-152, especially 
Theorem 8.9.4]. It also turns out that the Moore determinant is the same as the Jor- 
dan determinant on the Euclidean Jordan algebra Herm(n, H): see [561 Exercise II. 7 
(pp. 39-40), Exercise III. 1 (p. 58), and pp. 84 and 88]. This equivalence is what allows 
us to formulate Theorems 11.2( c) and II. 3( c) in terms of the Moore determinant, while 
the key Theorem 11.41 is stated and proven in terms of the Jordan determinant. 

Proposition A. 3 (Properties of the Moore determinant) 

The map Mdet: Herm(n, HI) — > C has the following properties: 

(a) Mdet is a homogeneous polynomial of degree n, with integer coefficients^ 

(b) Mdet / = 1. 

(c) Mdet is real-valued. 

(d) If M is a complex hermitian matrix, then its Moore determinant equals its 
ordinary determinant. 

(e) If X E M nxn (not necessarily hermitian), then Mdet(XX*) = det($(X)) = 
det(^(X)) > 0. 

(f) IfM e Herm(n,H) andX G H nxn ; then Mdet (X MX*) = Mdet(XX*) Mdet(M). 

(g) (restricted Cauchy— Binet formula for the Moore determinant) If X E 

Hjmxn with m <n, and D e IR nxn is a real diagonal matrix, then 

Mdet(X£>X*) = Mdet[X* 7 (X* 7 )*] det(D n ) . (A.14) 

JC[n] 
|/|=m 



37 More precisely, if for M = (rriij) € Herm(n,H) we write my = onj \ + fiifi + -fijj + Sijk with 
a%j , Pij , , Sij e R, then Mdet(Af) is a homogeneous polynomial with integer coefficients in the 
variables {au}i<i< n and {a^ ■, (3 i:j ■, Sij}i< i<: j< n . 
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(h) (restricted multiplicativity) If M,N <E Herm(n,H) with MN = NM, then 
Mdet(MX) = Mdet(M) Mdet(X). 

Proof, (a) is an immediate consequence of Lemma [A. 1( a) together with the defini- 
tions (IA~T3|) and (lA~5l) / fTA~8|) . 



(b) is likewise immediate from Lemma [A. 1( b). 

(c) Note that for any Y G $(Herm(n, H)) we have JYJ T = Y, hence 



pf(y) = pf(y) = pf(jyj i ) = (detj)pf(y) = P f(y) (A.15) 

by Lemma IA.l( d) 

mnlpv Viprmi1"inn mnfriv iV\pr\ <3}l A/f i = I 

-M 



(d) If M is a complex hermitian matrix, then <£>(M) = ( ^ ) and Mdet(M) 



pf (-m T) =detMb y ra>- 



(e,f) We have Mdet(XMX*) = pf[J#(XMX*)] = pf [J*(X)*(M)*(X*)] = 
pf [$(X) J$(M) J$(X*)] = pf [-$(X) J$(M) J$(X) T ] = pf[($(X) J) $(M) ($(X) J) T ] 
det($(X)J) pf($(M)) = det($(X)) Mdet(M). Applying this with M = I, we 
see that Mdet(XX*) = det($(X)) = det(#(X)). Therefore, Mdet(XMX*) = 
Mdet(XX*) Mdet(M). 

We will prove in Lemma lA.4( c) below that det(\I/(X)) > 0. 

(g) Define \l/ mn : H mxn — >• C 2mx2n by the same formula (IA.5|) : then it is a homo- 
morphism in the sense that if M G M mxn and X G H nxp , then 

^ mp (MX) = ^ mn (M) V np (N) . (A.16) 

We also have, for M G H mxn , 

^ nm (M*) = -J n ^ mn (M) T J m . (A.17) 

Then we have 

Mdet(XDX*) = pf[J m * mm (XL>X*)] 

= pf[J m * mn (X)* nn (L>)# nm (X*)] 

= pf[J m y mn (X)[-J n $(D)}i-J n y mn (X) T J m }} 

= pf[(J m y mn (X)J n ) $(£>) (J m V mn (X) J n ) T ] 

= Yl det[(J m ^ mn (X)J n ),/] (pf $(£>)//) • (A.18) 

/ C [2n] 
|J| = 2m 

Now write I = K U L with X C [n], L C [2n] \ [n] and |X| + \L\ = 2m. We claim 
that pf Q(D)n = unless L = K + n. To see this, note first that for any complex 
matrix D (diagonal or not) we have 

9(D) = J*(D) = °y (A.19) 
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If in addition D is diagonal and L ^ K + n, then &(D)n has a zero row (and a 
zero column), so that det = 0. If furthermore D is real, then Q(D)u is 

antisymmetric and (pf $(-D)//) 2 = det &(D)u = 0. So only the terms L = K + n 
survive in the sum over I. 
Then in this case we have 

pf = pf (_^ A . D q*) = det£ >*x ( A - 2 °) 

by (|A.12p . Moreover, a simple calculation shows that 

{Jm^ mn{X)J n )- k l = J m ^ mm^X+K) J m (A. 21) 

and writing Y = X+ K we see easily that det(Jtf (F) J) = det($(F) J) = det($(F)) = 
Mdet(FF*) by part (e), hence 

det[(4$ mn (X)J n ) iJ ] = Mdet[X„ x (X^)1 . (A.22) 

(h) If M,N G Herm(n,H) with MJV = NM, we have M G Herm(n,H), 
so its Moore determinant is well-defined; then [Mdet(iWTV)] 2 = [pf (f(MiV))] 2 = 
det($(MJV)) = det (MAO) = det(*(M)*(iV)) = det(*(M)) det(*(A0) = 
det($(M)) det($(iV)) = [pf($(M))] 2 [pf ($(A0)] 2 = [Mdet(M)] 2 [Mdet(A^)] 2 . Now 
apply this with M, N replaced by M A = (1 - X)I + AM, A" A = (1 - X)I + XN where 
A G R. Then P(A) = Mdet(M A A" A ), Q(\) = Mdet(M A ) and R{\) = Mdet(A" A ) are 
polynomials satisfying P(A) 2 = Q(A) 2 P(A) 2 and P(0) = Q(0)P(0). It follows that 
P(A) = Q(A) = R(X) = 1. Now evaluate at A = 1. □ 



Lemma A. 4 (a) Let M G C nxn . Then the eigenvalues of MM are real or come in 
complex- conjugate pairs; and the negative real eigenvalues have even algebraic 
multiplicity. 



(b) If M G C nxn and k>0, then det(«I + MM) > 0. 

(c) If A, Be C nxn , then 

-B A 



del ( \ % I > . (A.23) 



Proof, (a) It is well known that AB and BA have the same characteristic poly- 
nomial. Therefore, MM and MM have the same characteristic polynomial, i.e. the 
coefficients of this characteristic polynomial are real. It follows that the eigenvalues 
of MM are real or come in complex-conjugate pairs. 

If MM has n distinct eigenvalues, then we can argue as follows: Suppose that 
MMx = Xx with x ^ 0. Then by taking complex conjugates we have MMx = 
Ax; and left-multiplying by M we obtain MM(Mx) = X(Mx). If A is real, the 
distinct-eigenvalues hypothesis implies that Mx = fix for some // G C By complex- 
conjugating we get Mx = fix, and left-multiplying by M we obtain MMx = JiMx = 
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|/i| 2 x. Since x ^ it follows that A = \fi\ 2 > 0. Therefore all real eigenvalues of MM 
are nonnegative. 

Let us now observe that the matrices M for which MM has n distinct eigenvalues 
are dense in C nxn . To see this, let us consider the real and imaginary parts of the 
entries of M to be 2n 2 distinct indeterminates, and let us form the discriminant of 
the characteristic polynomial of MM. This is a polynomial (with real coefficients, 
though we don't need this fact) in the 2n 2 indeterminates, and its zero set in M. 2n 
corresponds precisely to the matrices M G <C nxn for which MM does not have n 
distinct eigenvalues. Now the zero set of a polynomial in R N is either all of M or 
else a proper subvariety (which in particular has empty interior). Since there do exist 
matrices M G C nxn for which MM has n distinct eigenvalues, it follows that such 
matrices form a dense open set in C nxn . 

Since the eigenvalues depend continuously on the matrix, it follows by density and 
continuity that, for arbitrary M G C nxn , the negative real eigenvalues of MM have 
even algebraic multiplicity. 

(b) follows immediately from (a), using the fact that the determinant of a matrix 
is the product of its eigenvalues taken with their algebraic multiplicity. 

(c) If A is invertible, then we have (writing C = A~ X B) 



(A B_\ (A 0_\ f I_ C\ [A \ f I C\ fl + CC \ 
\-B A J \ Aj \-C I J \0 Aj \0 I J \ -C I ) 



which is nonnegative by part (b). Since the invertible matrices are dense in C nxn , 
the general result follows by continuity. □ 

Remarks on Proposition IA.31 I. The assertion in Proposition IA.3( e) that 
Mdet(X X*) > can alternatively be proven as follows: Observe that XX* is positive- 
semidefmite and hence has nonnegative real eigenvalues Corollary 5.3 and Propo- 
sition 5.2]; so by the spectral theorem [571 Theorem 3.3 and Proposition 3.8] we can 
write XX* = UDU* with UU* = U*U = I and D diagonal with nonnegative real 
entries; it then follows from part (f) that Mdet(XX*) = Mdet(J7J7*) Mdet(D) = 
Mdet(D), which is > by part (d). 

2. The restriction in part (g) to real diagonal matrices is essential: indeed, already 
in the first nontrivial case m = 1, n — 2 one sees that 



(A.24) 



from which it follows that 




(A.25) 



rar + rqs + sqr + sbs 



(A.26) 



while the right side of flA.14j) is (rr)a + (ss)b. 
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To see why part (g) states the appropriate Cauchy-Binet formula for the Moore 
determinant, recall first the ordinary Cauchy-Binet formula for matrices over a com- 
mutative ring R: for X G R mxn , M G R nxn and Y G R nxm , we have 

det(X MY) = (det X^) (det M u ) (det Yj*) . (A.27) 

I,JC [n] 
\I\ = \J\ = m 

But in our case the determinant is defined only for hermitian matrices, so we need to 

(i) require Y = X* and M = M* so that XMY is hermitian, 

(ii) restrict to I = J so that Mjj is hermitian, and 

(iii) rewrite (det X^)(det(X*) Iir ) as det[X itI (X i<I )*). 

But restricting the sum to / = J is correct in the commutative case (e.g. if all the 
quaternionic matrix elements happen to belong to M) only if det Mjj = whenever 
I 7^ J, and this happens if and only if M is diagonal (the "only if" is seen by 
considering m = 1). So we really need to require that M be a real diagonal matrix, at 
least when m = 1. On the other hand, in the square case m = n we have I = J = [n] 
automatically, so this case does not require M to be diagonal [cf. part (f)]. 

The special case D = I of f lA.14j) was proven by Liebendorfer [961 Theorem 1]. 



We do not know whether the general result is new. 

3. Regarding part (h), note that the product MN is hermitian if and only if M 
and N commute, so this hypothesis is required for Mdet (MiV) to be well-defined. 

Part (h) can alternatively be proven by invoking the spectral theorem for com- 
muting hermitian quaternionic matrices [57] Theorem 3.3 and Propositions 3.5, 3.6 
and 3.8]: if M,N G Herm(n, M) with MN = NM, then there exist a unitary ma- 
trix U G H nxn (i.e. UU* = U*U = I) and real diagonal matrices D,E such that 
U*DU = M and U*EU = N. The result is then an easy consequence of parts (d) 
and (f). □ 

Remarks on Lemma IA.41 1. Lemma lA.4( a) goes back at least to Asano 
and Nakayama [8, Satze 19-21], who show that MM is similar to the square of a 
real matrix. It is part of more general theorem of Youla [T4"8] Theorem 1] that gives 
a normal form for matrices M G C nxn under unitary congruence, i.e. M = U^HU 
with U unitary and £ of a special form. See also [761 P- 147], [EH Sections 4.4 and 
4.6, especially pp. 252-253], and [58]. Lemma lA.4( c) goes back at least to Zhang 
|151[ Proposition 4.2] — with the same proof as given here — but is probably much 
older@ 

2. In Lemma IA.4( a). it is possible for MM to have negative real eigenvalues: 
consider, for instance, the 90° rotation matrix M = ( ^ J ^ . It is also possible for 

MM to be nondiagonalizable: consider, for instance, M = ( J ^ j with e G R\ {0}. 

38 Roger Horn informs us that the forthcoming second edition of [81] will contain a vastly expanded 
treatment of these topics. 
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3. The assertions in Lemma E3](a) about the eigenvalues of MM are best possible, 
in the sense that every n-tuple (Ai, . . . , A n ) G C n satisfying these conditions is the set 
of eigenvalues for some matrix MM (indeed, one that is real and diagonalizable). To 
see this, it suffices to consider lxl matrices M = (A) with A > 0, and 2x2 matrices 

(cos 9 sin 9 \ 
. „ „ with A > and 9 G R, and then form direct sums. 
— sin 9 cos 9 J 

4. If we consider the real and imaginary parts of the entries of M G C nxn to 
be 2n 2 distinct indeterminates, then Lemma [A. 4( b) asserts that det(J + MM) is a 
nonnegative polynomial on ¥L 2n . It would be nice to have a more direct proof of this 
fact. For instance, is this nonnegative polynomial actually a sum of squares? And 
likewise for Lemma [A. 4( c). □ 



We also have a restricted formula for the Moore determinant of the sum of two 
matrices: 

Proposition A. 5 [5, Proposition 1.1.11] If M G Herm(n,H) and D G R nxn is 
a real diagonal matrix, then 

Mdet(£> + M) = ( det £>//)(Mdet M/ c/c ) . (A.28) 

IC[n] 

Remark. For matrices over a commutative ring we have an unrestricted formula 
for det (A + B), which involves a double sum over subsets /, J (Lemma 16. ip . But 
in the quaternionic case we must restrict to I = J in order to ensure that Ajj and 
Bjcjc are hermitian; and to justify this restriction we must require that at least one of 
the matrices A, B be diagonal. Therefore only the restricted formula (Corollary 16. 2p 
generalizes to quaternions. □ 



For matrices over a noncommutative division ring (such as the quaternions), one 
can in principle define four ranks: the left and right row-ranks and the left and right 
column-ranks. In general the left row-rank equals the right column-rank, and the 
right row-rank equals the left column-rank, but these two numbers need not be equal: 

fl i\ 

for instance, the matrix [ . ^ J has left row-rank and right column-rank 2, but right 

row-rank and left column-rank 1. Because matrix multiplication always forms right 
linear combinations of columns and left linear combinations of rows, it is left row-rank 
= right column-rank that is well-behaved with respect to multiplication; we therefore 
follow Jacobson [S71 pp. 22, 51] in focussing attention on these two, which we call the 
"rank" tout court. 

The rank of a general quaternionic matrix can be characterized as follows: 

Proposition A.6 [1171 Theorem 8] fT5T], Theorem 7.3] M G HP xn has rank r 
if and only iffy(M) [or equivalently $(M)/ has rank 2r. 
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For hermitian quaternionic matrices we conjecture an alternate characterization in 
terms of principal minors, which if true would extend a well-known result for complex 
hermitian matrices^!: 

Conjecture A. 7 Let M G Herm(n, H) have rank r. Then: 

(a) All principal minors of M of order r (i.e. the quantities MdetM/j with \I\ —r) 
have the same sign, and at least one of them is nonzero. 

(b) All principal minors of M of order s > r vanish. 

We also have an analogue of Proposition 16.31 for hermitian quaternionic matrices: 
Proposition A. 8 Let A\, . . . , A^ G Herm(n, M), and define 

P( Xl ,...,x k ) = Mdet (y^XiAA (A.29) 

for xi, . . . , Xk G R. Then P is a homogeneous polynomial of degree n with real coeffi- 
cients. Furthermore, the degree of P in the variable Xi is < rank(Aj). 



Proof. It follows immediately from the definition of Mdet that P is a homogeneous 
polynomial of degree n that is real- valued when real. Moreover we have 

P(xi,...,x fc ) 2 = det$^x^j = det^^$(A)J • (A.30) 

It then follows from Proposition 16.31 that the degree of P 2 in the variable Xi is at most 
rank($(Aj)), which equals twice rank(v4j) by Proposition IA. 61 □ 
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